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Abstract 

This report contains both a review of recent approaches to supersymmetric lattice 
field theories and some new results on the deconstruction approach. The essential 
reason for the complex phase problem of the fermion determinant is shown to be 
derivative interactions that are not present in the continuum. These irrelevant opera- 
tors violate the self-conjugacy of the fermion action that is present in the continuum. 
It is explained why this complex phase problem does not disappear in the continuum 
limit. The fermion determinant suppression of various branches of the classical mod- 
uli space is explored, and found to be supportive of previous claims regarding the 
continuum limit. 



*gied.t®physics .umn. edu 



Contents 



1 Introduction 



1.1 Motivations 



1.2 Tlie difficulty 



1.3 Summary 



Various approaches to supersymmetry on the lattice 



2.1 Deconstruction models 



2.2 Q-exact quantum mechanics 



2.3 Q-exact 2d Wess-Zumino models 



2.4 Supersymmetric nonlinear sigma model; 



2.5 Twisted sreo metrical SYM 



2.6 Q-exact compact super- Yang-MilL 



2.7 Perfect action model 



I 



2.8 Super-renormalizable lattice field theories 



2.9 Nonlinear O modelsl 



2.10 Chiral fermion models of 4d TVT = 1 SYM 



2.11 Chiral fermion formulations of the 4d Wess-Zuminomodel 



2.12 Twisted superspac^ modelJ 

2.13 Brute force model; 



2.14 Transverse lattice models 



2.15 Hamiltonian approaches 



Deconstructed (2.2) 2d SYM 



3.1 The continuum theory 



3.2 Construction of the daughter theory 



3.3 Classical approach to the continuum theory 



3.4 Stabilization of the a-configurationi 

3.5 Renormalization] 



3.6 The fermion determinant 



3.7 Comparison to U(k) and SU(k) matrix model 



3.8 Source of the complex phase problem 



4 Deconstructed lattice (4.4) 2d SYM 



4.1 The continuum theory 



1 



4.2 The mother and daughter theorie! 



4.4 The semi-classical analvsii 



4.4.1 Undeformed theory 



36 



4.3 Construction of the lattice theorvl 38 



39 
39 



4.4.2 



Deformed theorvl 42 



4.4.3 



Inhomogeneous minima of the deformed theorvl 44 

4.4.4 Determinant on moduli spacel 46 



4.5 Simulation analvsis of the dvnamic lattice spacing! 48 

4.5.1 Scaling! 49 



4.5.2 Problems with the complex phase| ■ ■ ■ ■ 49 



4.5.3 Results for the phase-auenched ensemble! 51 

4.5.4 Phase reweightind 52 



5 Conclusions 



A A toy (2.2) theory 



B Boundedness of the partition function 



C Deformed fermion matrix 



53 
55 
56 
57 



D Conjugation of the action and properties of the fermion determinant! 57 



D.l An example from kindergarten 


D.2 n-component complex Grassmann 




D.3 Imaginarv-time continuation of hermiticitv 





59 
60 



E AT = 4 moduli space 



60 



1 Introduction 

Several motivations exist for efforts to formulate supersymmetric field theories on a 
lattice. These will be summarized in ^1.11 It is difficult to formulate these theories 
in such a way as to avoid fine-tuning of counterterms. This is briefly discussed in 
§1.21 The points that I discuss in §l.ll - §1.2l are all well-known; experts may prefer to 
skip directly to §1.31 There I briefly summarize the content of this report, with an 
emphasis on new results that have not appeared elsewhere. 



2 



1.1 Motivations 



One is just to furnish a consistent nonperturbative definition for a given supersym- 
metric field theory. A quantum theory that is well-defined in perturbation theory 
may nevertheless have a nonperturbative anomaly. A classic example is, for example, 
the Witten anomaly [1]. A more recent, exotic example is [2]. With respect to func- 
tional integrals that appear in quantum field theory, the lattice approach provides 
a rigorous, nonperturbative way to give them concrete meaning, and to study their 
behavior. If there is a nonperturbative anomaly associated with one of the pertur- 
bative symmetries of the continuum theory, it would be interesting to see how that 
surfaces in the lattice theory. It has even been argued that it is possible to have a 
supersymmetry anomaly [3], which would be particularly interesting to study in the 
lattice formulation of such theories. 

A second motivation is that some nonperturbative aspects of supersymmetric field 
theories are not calculable (in a quantitative sense) by the usual techniques. One 
prime example occurs in cases involving chiral superfields. So-called F-term parts of 
the action {holomorphic data) are necessarily obtained from a superpotential, which 
is a holomorphic function of the chiral superfields.^ This part of the action is typ- 
ically subject to powerful nonrenormalization theorems. Together with symmetry 
constraints, this often allows one to draw conchisions about the effective superpoten- 
tial; the validity of these conclusions sometimes extends beyond perturbation theory. 
However, parts of the action that have a nonholomorphic origin, so-called D-term 
contributions {nonholomorphic data), can undergo significant renormalization when 
running couplings turn strong. Quite often, the continuum methods that are available 
(e.g., instanton calculus and R'^ x formulations) are not able to make definitive 
statements about this renormalization. As a consequence, quantitative or even quali- 
tative features of the low energy effective theory are not under theoretical control, to 
the extent that they depend on nonholomorphic data. In phenomenological applica- 
tions, this lack of control over nonholomorphic data is distressing, since the so-called 
supersymmetry- breaking soft-terms that partly determine supcrpartner spectra and 
couplings will depend on parts of the effective action that involve nonholomorphic 
quantities, such as the Kdhler potential. See for example [5,6]. Lattice Monte Carlo 

^In this discussion, the action may include source terms for elementary fields or composite oper- 
ators. In this case, holomorphic data may include composites than can be viewed as source terms in 
the superpotential. Equivalently, a superspace description of Green functions could be used; Green 
functions that are F-terms then comprise holomorphic data. See for example Chs. 9-10 of [4]. 
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simulations would, as a first step, give us a handle on the spectrum of light states. 
Through the introduction of source fields, it might even be possible to constrain 
the leading terms in the effective Kahler potential, when the lattice data is com- 
bined with what is known about the effective superpotential by continuum methods. 
This is because both the Kahler potential K and superpotential W play a role in 
the scalar potential: V = K^^WkWi, where K^^ is the inverse of the Kahler met- 
ric Kj^i = d'^K/d(j)^d(f)^ and Wk — dW/dcpk, — {Wi)*. Thus, lattice simulations 
have the potential to teach us something about nonperturbative renormalization of 
nonholomorphic quantities. 

A third motivation, which is closely related to the second, is that we would like 
to improve our understanding of dynamical supersymmctry breaking. 4d supersym- 
metric gauge theories play an important role in phcnomcnological applications. A 
strongly coupled sector is often invoked in the mechanism that splits the superpart- 
ners from the observed Standard Model spectrum. As examples, see [7] and references 
therein for a review of supergravity mediation of gaugino condensation^ and see [8] for 
a review of the gauge mediation scenario. Any improvement of our understanding 
of the strongly coupled gauge theories that arc invoked in these proposals would be 
helpful. If a good formulation of super- Yang- Mills (SYM) on the lattice is discov- 
ered, it is likely that interesting, well-posed questions could be studied through lattice 
simulations. 

A fourth motivation is that for certain super symmetric theories, quite a lot is 
known precisely, and it would be instructive — for the development of lattice theory 
and techniques — to reproduce these results on the lattice, in simulations or by analytic 
methods. For example, it would be very nice to reproduce the exact mass formulae, 
confinement and chiral symmetry breaking of pure J\f — 2 SYM that was derived by 
Seiberg and Witten in [9, 10]. 

A fifth motivation has to do with the evolving understanding of the relation- 
ship between SYM and theories of quantum gravity; in particular, string/M-theory. 
A nonperturbative formulation of superstring theory in general backgrounds is still 
lacking, in spite of the recent successes of M(atrix) theory [11,12], the AdS/CFT 
correspondence [13-15], the PP-wave hmit [16-19], etc. However, through the studies 
just cited, nonperturbative formulations in special backgrounds have emerged, and it 
is of considerable interest to study them in relation to SYM on the lattice for various 
numbers of spacetime dimensions. In particular, both the Matrix theory formula- 
tions of string/M-theory, and the AdS/CFT correspondence, are expressed in terms 
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of quantum theories of dimensionally reduced SYM. The nontrivial vacuum dynamics 
takes on a gravitational meaning. 

1.2 The difficulty 

The difficulty in formulating supersymmetric field theories is well-known, and was 
pointed out long ago by Dondi and Nicolai [20]. The supersymmetry algebra fits into 
the larger, super- Poincare algebra. In particular, the supersymmetry algebra closes on 
the generator of infinitesmal spacetime translation. Since the spacetime translation 
group is explicitly broken by the discretization, the supersymmetry algebra invariably 
must be modified on the lattice. The obvious option is to have it close on discrete 
translations. However, the Leibnitz rule does not hold on the lattice. A result of 
this is that even if the supersymmetry algebra holds on linear expressions of lattice 
fields, it is violated for more general polynomials. It then follows that interacting 
supersymmetric actions will not be invariant w.r.t. the lattice supersymmetry. 

Since the violation of the Leibnitz rule is an 0{a) effect, where a is the lattice 
spacing, the non-invariance of the lattice action w.r.t. lattice supersymmetry is like- 
wise 0{a). Classically this disappears in the continuum limit. Unfortunately, in the 
quantum theory the 0{a) violations can be overcome by UV divergences. Essentially, 
any non- irrelevant operator that is allowed by the symmetries of the lattice action 
will be radiatively generated. The lack of exact supersymmetry in the lattice theory 
generically leads to supersymmetry violating operators in the renormalized theory. 
These can be subtracted off with counterterms. However, this calls for a fine-tuning 
of the counterterm coefficients in order to achieve the desired continuum limit. A non- 
peturbative method of fine-tuning is required if one wishes to go beyond perturbation 
theory. 

Symmetries of the lattice theory can sometimes overcome this difficulty. This 
approach is discussed in much more detail below. In particular, see ^ below. 

1.3 Summary 

One method of building lattice models with some exact supersymmetry is based 
on dimensional deconstruction [21,22]. In I will focus on this approach as it 

applies to 2d SYM with extended supersymmetry; furthermore, I will concentrate on 
the formulation that leads to a Euclidean action formulation of the lattice theory, as 
opposed to a Hamiltonian formulation, where time is left as a continuous variable. 
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The 2d models are convenient in that they provide a simple context to illustrate the 
methods involved. Furthermore, these 2d SYM models are relevant to the AdS3/CFT2 
correspondence — an extension of the original AdS5/CFT4 ideas. (Here, subscripts 
indicate spacetime dimensions.) Generalizations of the deconstruction approach to 
higher dimensions will only be discussed parenthetically, since the techniques are the 
same as in the 2d case. 

The 2d examples already illustrate the principal point: sufficient cleverness in 
constucting the lattice formulation can overcome the apparent difficulties in obtaining 
the desired continuum limit. Unfortunately, UV behavior gets worse for d > 2. This 
tends to make the avoidance of unwanted renormalizations — i.e., those that produce 
supersymmetry breaking non-irrelevant operators — more challenging, especially in 
the case of most interest, d = 4. For d > 4, a lattice SYM formulation would only 
be useful as an effective theory that captures nonperturbative effects with a scale far 
below the cutoff, since these theories are nonrenormalizable.^ Due to infrared freedom 
in these theories, continuum methods are already reliable in this regime — although 
lattice analysis could nevertheless prove instructive, d > 4 theories probably require 
something more along the lines of string/M-theory in order to make sense as anything 
other than a low-momentum effective theory. 

The deconstructed spacetime lattice constructions are all arrived at by orbifold 
projections of SYM matrix models; i.e., in each case we quotient a Od matrix model 
(SYM with all spacetime dimensions reduced out) by some discrete symmetry group 
of the theory. Degrees of freedom that are not invariant with respect to the orbifold 
group are projected out. This approach was introduced for a Euclidean spacetime 
lattice formulation^ by Cohen et al., whom I refer to hereafter as CKKU [24]. In this 
report I will examine both the (2,2) and (4,4) 2d SYM constructions.^ 

I now summarize the remainder of this report: 

• Before delving into the details of the CKKU models, I first provide in ^ a 
summary of various approaches to supersymmetry on the lattice, both non- 
gauge and gauge models. 

• In Pl discuss aspects of the (2,2) 2d SYM construction of CKKU [24]. In 

^In this regard, the situation would be analogous to lattice heavy quark effective theory or lattice 
nonrelativistic QCD. 

•^Kaplan et al. had previously formulated a spatial lattice by a similar method [23]. 
^Here and below, (m,n) refers to m left-chirality and n right-chirality supercurrents (equivalently, 
m left-moving and n right- moving supercurrents). See Appendix 1X1 
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particular, I describe how the phase of the fermion determinant is generically 
complex, in a boson configuration dependent way [25]. In this report I present a 
new, more fundamental understanding of this problem. The origin of this phase 
is shown to be 0{a) irrelevant, lattice derivative operators that violate the self- 
conjugacy that is enjoyed by the continuum fermion action. Due to the presence 
of modes with 0{a~^) momenta in the fermion determinant, these derivative 
operators make a non-negligible contribution, even in the a —>■ limit. 

• In gll delve into the (4,4) 2d SYM construction of CKKU [26]. In addition to 
other matters, a discussion is given of the classical moduli space of the lattice 
theory, and the effects of the magnitude of the fermion determinant for fiuctua- 
tions about points on this moduli space. Although some part of this discussion 
has appeared elsewhere [27], a more thorough analysis is contained here — with 
results that were not previously noted. Some simulation results in the phase- 
quenched ensemble are briefiy described, including my attempts to reweight by 
the phase of the fermion determinant — which turned out to be unsuccessful. It 
is shown, however, that phase- quenched averages support the general picture of 
how the continuum limit emerges. 

• A brief set of conclusions is contained in I also summarize related research 
that I now have in progress. 

• App. ^ describes the simplest type of 2d field theory with (2,2) supersymme- 
try. This is provided as an illustration of the basic structure of this sort of 
supersymmetry, and the logic of the notation "(2,2)." 

• App. El briefiy reviews what is known about the boundedness of the matrix 
partition functions that are used in the deconstruction approach to lattice SYM. 
It is described how noncompact fiat directions of the classical action are rendered 
effectively compact in the "quantum" theory — i.e., when integration over modes 
orthogonal to the moduli space is performed. 

• A result that was lacking in earlier works is presented in App. O Here I rigor- 
ously prove that the deformed fermion matrix correctly factors out ever-present 
fermion zeromodes that occur in the deconstructed models, and gives a projec- 
tion onto the orthogonal complement. 
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• In App. ini a review of generalized conjugation of fermionic theories is pro- 
vided. It is described how unitary conjugation maps may be exploited to prove 
important properties of the fermion determinant, such as reality. 

• App.|Elgives a brief proof of the description of the A/" = 4 4d SYM moduli space, 
since the same algebraic constraints must be satisfied in (4,4) deconstructed 
model that is described in ^of the main text. 

2 Various approaches to supersymmetry on the 
lattice 

Having motivated the study of lattice supersymmetry, I next give a brief survey of 
some of the more promising approaches that have appeared in the literature of late. 

We know from experience with other continuum symmetries (e.g., Euclidean ro- 
tation invariance and chiral symmetry) that exact lattice symmetries are often the 
key to obtaining the desired continuum limit. This is because the exact lattice sym- 
metries forbid unwanted operators: those that are non-irrelevant (according to naive 
power-counting) and would otherwise destroy the desired symmetry in the quantum 
continuum limit, absent fine-tuning of counterterms. Thus, exact lattice symmetries 
can in some cases guarantee that the long-distance Wilson effective action will have 
the desired continuum symmetries. 

Indeed, in some cases it has been shown how to construct lattice actions whose 
symmetry group forbids unwanted supersymmetry-breaking operators. In most cases, 
this occurs due to an exact lattice version of some part of the continuum supersym- 
metry algebra, in conjunction with other lattice symmetries. In some cases, there 
is only a partial success, in that the number of counterterms that must be adjusted 
is reduced but not zero. In that case, some tuning of bare parameters or additional 
counterterms will still be required to achieve the desired continuum limit. 

One should not lose sight of the fact that a good continuum limit typically arises 
from the combined constraints of all lattice symmetries — the space group of the lat- 
tice, discrete internal symmetries, a subalgebra of the supersymmetry, gauge symme- 
try, etc. I would also like to point out that forbidding unwanted operators may or 
may not be sufficient to guarantee the correct continuum limit. An example from 
4d phenomenology is the following. Generic 2-Higgs Doublet models and the Min- 
imal Supersymmetric Standard Model differ, among other things, in one important 
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respect: in the former case the strength of the quartic couphng in the Higgs potential 
is arbitrary, whereas in the latter case it is related to the gauge couplings in a precise 
manner.^ More generally, supersymmetric theories are typically more constrained 
than the most general theory containing the same operators and symmetries other 
than supersymmetry. These additional constraints would need to be guaranteed by 
the exact lattice symmetries if one is to be entirely successful in avoiding fine-tuning 
of counterterms. 

I now briefly summarize the existing threads of construction, based mostly on 
exact lattice symmetries. 

2.1 Deconstruct ion models 

These are lattice formulations of super- Yang-Mills that follow an approach based on 
dimensional deconstruction [21,22]. In deconstruction, dimensions are replaced by 
product groups <S>me{sites}'^rn connected by representations^ TZm,i = {Rm, Rm+v(.i)) 
that link the factors Gm, just as in lattice gauge theory; however, in deconstruction 
the lattice spacing that defines the continuum limit is determined by a background 
value of the link fields.^ Examples in the literature include: a spatial lattice [23]; 
Euchdean spacetime lattices [24,26,30]; and, lattices of just one spatial dimension 
[31,32] . Note that I have not listed the constructions that were the original motivation 
for dimensional deconstruction: theories with spacetime dimension d > A. That is 
because such theories are nonrenormalizable and do not have a continuum limit, in 
the lattice field theory sense. Here our interest is only in theories where the lattice 
cutoff can ultimately be removed. The deconstruction approach is the subject of the 
bulk of this report. Fine-tuning is very much reduced in these constructions, and in 
some cases nonexistent. 

2.2 Q-exact quantum mechanics 

In this case. Id supersymmetric quantum mechanics [33,34] is obtained from a Q-exact 
action (i.e., S = QX). Here, Q is an exact lattice supercharge that is nilpotent, — 

^Here I refer to the quartic operators that are allowed by supersymmetry, and that arise from 
D-terms. 

^Herc, {u*^*)} is a set of constant vectors, composed of elementary lattice vectors. 

^It is amusing to note that product group theories with this sort of spontaneous symmetry 
breakdown to a diagonal subgroup were considered a long time ago, but without the dimensional 
deconstruction interpretation [28,29]. I thank Marty Halpern for bringing this to my attention. 
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0, once auxiliary fields are introduced. In derivatives are realized through discrete 
difference operators; with respect to a discrete approximation of the continuum theory 
superalgebra, = is a subalgebra. It is trivial that the action is "supersymmetric" 
(i.e., invariant w.r.t. Q): QS — Q'^X — 0. Models in this class have been derived 
from: a Nicolai map [33,34], a relationship to topological field theory [35,36], and 
superfield techniques [37]. Fine-tuning of counterterms is absent in the continuum 
limit. This has been shown in perturbation theory and at a nonperturbative level [37]. 
By contrast, in a naive discretization a counterterm must be introduced [38]. 

2.3 Q-exact 2d Wess-Zumino models 

In this case the lattice theory is supposed to provide a Landau-Ginsburg description 
of the minimal discrete series oiAf — 2 superconformal field theories. The continuum 
Landau-Ginsburg effective theory is nothing but the 4d — > 2d dimensional reduction 
of 4d Wess-Zumino models [39] . 

These Q-exact lattice actions were first formulated in [40, 41] using Nicolai map 
[42] methods, relying on earlier Hamiltonian [43] and continuum [44] studies that 
also utilized the Nicolai map. This Nicolai map approach was also described in 
[34,45]. The same class of lattice actions was obtained by superfield techniques 
in [37] . Once auxiliary fields are introduced, the lattice action takes a Q-exact form, 
as was emphasized in topological interpretation of [35, 36] and the lattice superfield 
approach of [37] . Here again, there is a nilpotent subalgebra Q"^ — that is preserved 
in the formulation with auxiliary fields. 

Detailed studies of the spacetime lattice system have been performed: in [46] by 
stochastic quantization methods; in [34,47] by the Monte Carlo simulation approach; 
and in [37] by lattice perturbation theory. 

It has been shown in [37] that the continuum limit of the lattice perturbation 
series is identical to that of the continuum theory, due to cancellations that follow 
from = 0. Thus, the Q-cxact spacetime lattice has behavior that is similar to 
what was found on the Q, <5^-preserving spatial lattice in [40]. However, it was also 
shown in [37] that the most general continuum effective action that is consistent with 
the symmetries of the bare lattice action is not the (2,2) 2d Wess-Zumino model. 
Thus the nonperturbative studies [34,46,47] cited above are important to settling the 
question as to whether or not the correct continuum limit is obtained. 

The spectrum degeneracy and approximate satisfaction of the supersymmetry 
Ward-Takahashi (WT) identities that was observed in the studies [34,46] indicate that 
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the correct continuum limit is obtained without fine-tuning, even at a nonperturbative 
level. In [47] the R-symmetries of this lattice system were studied. It was found that 
the continuum symmetries do seem to be recovered in the nonperturbative, quantum 
continuum limit. In summary, all indications are that the lattice system has the 
correct continuum limit without any fine-tuning. 

2.4 Supersymmetric nonlinear sigma models 

Supersymmetric 2d nonlinear sigma models have been attempted on the lattice. Here 
both topological [35,36] and superfield [37] approaches have been applied. One finds 
that it is necessary to give up either the exact nonlinear symmetry or the exact lattice 
supcrsymmetry, assuming the nonlinear symmetry group is nonabelian. It remains 
to be seen which is the better option and whether or not fine-tuning of counterterms 
is required. The nonlinear symmetries of the 2d sigma model are important to the 
renormalizability of the theory, so sacrificing them could create severe problems with 
the continuum limit. Encouraging analytic results and some Monte Carlo simulations 
of the CP^ model have been reported in [48,49]. 

2.5 Twisted geometrical SYM 

These are formulations that rely on the geometrical picture that relates twisting 
of spinors and the Kahler-Dirac formalism. With the scalar supercharge Q that is 
exposed in this procedure, and all elementary degrees of freedom encoded in Kahler- 
Dirac fields, Catterall has been able to construct Q-cxact actions for super- Yang-Mills 
(SYM) with extended supcrsymmetry. The ideas behind Catterall's construction go 
back to his topological approach for non-gauge theories, described in [35,36]. This 
is because twisting is an integral part of the formulation of topological field theories. 
In [50] the jV = 2 4d SYM was formulated, whereas in [51] the = 4 4d SYM was 
constructed. 

2.6 Q-exact compact super- Yang-Mills 

These are the models of Sugino [52-55], where his approach is to modify the super- 
symmetry transformation in such a way that it has a simple action on the link fields 
that appear in a compact formulation of SYM. He achieves this by supplementing 
the lattice supcrsymmetry tranformation with 0( a) -suppressed operators. He then 
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defines the action as a Q-exact quantity, where = up to a gauge transformation. 
By this approach he is able to construct a variety of extended SYM theories with 
an exact lattice supersymmetry. He demonstrates the absence of fine-tuning at the 
perturbative level, in way that is quite similar to what was done in the deconstruc- 
tion models. Recently, there appeared a brief note that questioned Sugino's claim 
that he has been able to obtain a unique classical vacuum about which to define the 
continuum limit [56]. However, this note is wrong in its conclusions, due to several 
errors that are so elementary that I need not detail them here. 

My own studies have found that the fermion action is not hermitian, due to 
irrelevant operators, indicating a complex fermion determinant. I will present these 
results in a forthcoming article. 

2.7 Perfect action models 

Direct constructions in the spirit of the Ginsparg- Wilson relation, as suggested by 
Liischer [57], have been considered in a few articles [58-60]. However, the only cases 
where this has been successfully applied are theories with quadratic actions; i.e., free 
theories. Quadratic theories can be successfully latticized in a naive approach, if 
a judicious choice of discrete difference operators is made. Furthermore, a lattice 
approach is not needed for quadratic theories, since they can be solved exactly in the 
continuum. Thus, the perfect action approach has not yet proven useful for interesting 
(i.e., interacting) supersymmetric field theories on the lattice. 

2.8 Super-renormalizable lattice field theories 

In [61], Golterman and Petcher studied the A/" = 1 2d Wess-Zumino model. They 
showed by Reisz lattice power-counting [62] that the lattice perturbation theory is 
super-renormalizable, with all diagrams of UV degree D > appearing at 1-loop. 
Thus, the counterterms that were needed to guarantee the perturbative continuum 
limit could all be exactly determined by a straightforward 1-loop computation. 

In the proof that the supersymmetry WT identities were respected in the con- 
tinuum limit, Golterman and Petcher exploited a modified lattice supersymmetry 
algebra that contained a series of higher dimensional operators, suppressed by posi- 
tive powers of the lattice spacing. This sort of nonlocal, exact lattice supersymmetry 
was inspired by the earlier work of Banks and Windey [63]. The Golterman and 
Petcher method of analysis inspired the nonlinear Q models that will be discussed in 
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^ below. 

More generally, consider the case of any lattice field theory that is super-renormaliz- 
able according to the Reisz power-counting rules. It follows that such a theory allows 
for a counterterm adjustment that will guarantee the correct perturbative continuum 
limit. Exactly this sort of approach has been carried out for 3d A/" = 2 SYM [64]. 
Another interesting result is [65], where (2,2) 2d SYM is formulated with chiral lat- 
tice fermions (cf. ^2.1(J|) . so that the only counterterms are scalar masses. These are 
entirely determined at 1-loop since the theory is super-renormalizable. 

However, the crucial role of the lattice is supposed to be as a nonperturbative 
definition of the quantum field theory. One must ask: To what extent does good 
perturbative behavior lead to good nonperturbative behavior? Here, there is only 
anecdotal data from a few nonperturbative analyses of supersymmetric lattice theo- 
ries; these have been referred to above or below. It does seem that good pertubative 
behavior correlates with good nonperturbative behavior for these few simple models; 
for example, adjustment of counterterms as determined from a super-renormalizable 
lattice perturbation theory has been found to give good nonperturbative results in 
the naive latticization of supersymmetric quantum mechanics [38] . Also, some models 
where it was shown that no counterterms are required in perturbation theory have 
been studied by nonperturbative means, with the result that the correct continuum 
limit appears to emerge at the nonperturbative level as well: the examples here are 
the Monte Carlo studies of Q-exact lattice supersymmetric quantum mechanics [38] 
(cf. ^2.2j) . and nonperturbative studies of the Q-exact A/" = 2 2d Wess-Zumino model 
by stochastic [46] and Monte Carlo approaches [34,45,47] (cf. §2.3|) . But without a 
rigorous proof, one cannot say in what cases these types of positive results will per- 
sist. The best one can do is to test each theory nonperturbatively; say, by checking 
supersymmetric WT identities. 

2.9 Nonlinear Q models 

The work of Goltermann and Petcher [61] has recently inspired Bonini and Feo to 
formulate a nonlocal, exact lattice supersymmetry for a discretized version of the 4d 
Wess-Zumino action [66-68]. This action makes use of Ginsparg- Wilson fermions, 
so that the U{l)fi symmetry of the massless theory is respected without fine-tuning 
[69-71]. While the nonlocal exact lattice supersymmetry proposed by Bonini and 
Feo is a symmetry of the partition function, order by order, it remains to be seen 
what it has to do with symmetries of the tranfer matrix, and hence Hamiltonian. 
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Note that an operator that is nonlocal in the time direction touches all time-slices. 
As a consequence, if this operator generates a symmetry of the partition function, 
it does not necessarily imply that the operator commutes with the transfer matrix. 
If the operator could be shown to localize in the continuum limit, then it would 
indeed commute with the coarse-grained Hamiltonian in this limit. I am currently 
investigating these matters. 

The analysis of the continuum limit of perturbative lattice WT identities that 
was reported by Bonini and Feo in [68] indicates that the exact nonlocal lattice 
symmetry is related to a symmetry of the continuum Hamiltonian through 0{g^). 
However, according to Kikukawa and Susuki [72], Ref. [69] shows that there is already 
supersymmetry breaking at 1-loop order in the wavefunction renormalization. It is 
unclear to me how this can be reconciled with the findings of Bonini and Feo regarding 
the WT identities. 

2.10 Chiral fermion models of 4d A/" = 1 SYM 

Here a fair amount of work has been done. It began with the observation of Curci 
and Veneziano that for 4d pure M — 1 SYM, the violation of the supersymmetry 
WT identity is proportional to the mass- violation of the axial U{1)a anomalous WT 
identity [73]. This is easy to understand: if one assumes gauge and Lorentz invariance, 
the only non-irrelevant operator that can break supersymmetry for 4d = 1 SYM 
is a gluino mass.^ 

Thus the supersymmetric limit coincides with the chiral limit of the lattice theory; 
that is, the limit of vanishing renormalized mass. Due to the axial anomaly, in the 
case of SU{N) SYM, only a Z2N subgroup of U{1)a is a symmetry of the partition 
function. This is spontaneously broken to Z2 due to condensation of the gluino 
bilinear, (AA) 7^ in a 2-component Weyl spinor notation. Consequently, one can 
either (1) fine-tune the bare mass for a lattice fermion that explicitly breaks the 
Z2N chiral symmetry, or (2) use a chiral fermion formulation — where there is only 
multiplicative renormalization, and setting the bare mass to zero will correspond to 
the chiral limit. 

In case (1), Wilson fermions have been explored extensively (see for example 
[74-76] and references therein). The chiral limit corresponds to a first order phase 
transition that occurs as one varies the bare mass such that one passes through the 
^The gluino is the fermionic superpartner to the gluon. 
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critical mass where N vacua coexist. See the review [75] for an extensive discussion 
of this matter and recent findings in numerical simulations. 

In case (2), several formulations have been considered. Staggered (or, Kogut- 
Susskind) fermions [77-79], which possess symmetries that forbid additive mass renor- 
malization, were studied in [80]. Because of a 4-fold increase in modes, and because 
SYM has Majorana fermions, leading to a Pfaffian,^ an eighth-root trick was ap- 
plied to the staggered fermion determinant. This technique is controversial, so other 
formulations are of interest, as a cross-check. 

Of particular interest are chiral formulations that are collectively referred to as 
Ginsparg- Wilson fermions, since they satisfy the Ginsparg- Wilson relation [81], a 
lattice variant of chiral symmetry [57]. Domain wall fermions [82] yield a chiral for- 
mulation in the limit of an infinite 5th dimension. Monte Carlo simulations of A/" = 1 
4d SYM have been performed for domain wall fermions by Fleming et al. [83]. The 
cost of such simulations is high, so only limited results have been obtained. The 
gluino condensate was studied; however, the data that was obtained was not suffi- 
cient to extrapolate to the continuum limit. The overlap determinant gives an exact 
representation of the Weyl determinant in the continuum limit [84,85]. Since in 4d a 
Weyl fermion is equivalent to a massless Majorana fermion, the overlap determinant 
is guaranteed to yield the correct continuum theory. A lattice Dirac fermion that 
enjoys exact chiral symmetry was later obtained by Neuberger, based on the overlap 
Weyl determinant [86,87], and this is what is commonly referred to as an overlap 
fermion; it is constructed from the Wilson-Dirac operator. In the case of the overlap 
formulation, only analytical studies have been performed [88], due to the fact that 
simulations are extraordinarily expensive; significant supercomputer resources would 
be required to obtain meaningful results for the 4d A/" = 1 SYM theory. 

2.11 Chiral fermion formulations of the 4d Wess-Zumino 
model 

Some of the efforts in this direction, using Ginsparg- Wilson fermions, were mentioned 

in §2.91 above. Such models were first described by Fujikawa and Ishibashi in [69-71]. 

^Up to a 2-fold branch ambiguity, the PfafRan of a matrix is the square root of its determinant. 
The branch may be determined from the precise definition of the Pfaffian, which is well-known and 
I will not repeat here. In the cases where the Pfaffian does not have constant complex phase (or 
sign, if it is real) , the determination of and sample weighting by this factor often inhibits numerical 
studies. 
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The analysis of that construction by Bonini and Feo was aheady summarized in 
§2.91 A modified formulation was given by Kikukawa and Susuki [72]. In their 
constuction they have shown that the lattice chiral symmetry reduces the number of 
counterterms that are required, but does not eliminate them. In order to obtain the 
correct continuum limit, a tuning must be performed at each order in perturbation 
theory. In light of this, the seemingly contradictory claims of Bonini and Feo regarding 
the Fujikawa and Ishibashi construction need to be re-examined. 

2.12 Twisted superspace models 

In [89,90], models have been introduced in which superspace is given a noncommu- 
tative, lattice-geometrical interpretation. For example, the supercharges correspond 
to a sort of "square-root" of lattice translation. Key in this construction is the 
Kahler-Dirac formalism and complexification of the gauge group. Studies of dynam- 
ical properties of these lattice systems have not yet been performed. 

2.13 Brute force models 

Here one just identifies the most general renormalizable lattice action that contains 
the supersymmetric target, and fine-tunes the bare parameters to obtain the desired 
theory. An example is the work of Montvay for 4d A/" = 2 SU(2) SYM [91]. There the 
most general theory consistent with gauge invariance and parity (easily preserved on 
the lattice) is an SU{2) adjoint Higgs- Yukawa theory, which he formulates with Wilson 
fermions and compact gauge fields (link variables). He then explores perturbative RG 
fiows and identifies the trajectory that gives the desired theory. A nonperturbative 
extension of this work would be very interesting, due to the existence of exact results 
in the continuum [10]. 

2.14 Transverse lattice models 

In this class of models, lattice techniques are combined with light-cone quantization. 
The dimensions transverse to the light-cone are formulated on a lattice. This allows 
for more symmetry to be preserved (associated with the light-cone), and reduces the 
dimensionality of the lattice system that must be simulated. An example of this type 
of construction is [92]. 
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2.15 Hamiltonian approaches 

In this case, time is treated continuously, whereas spatial dimensions are discretized. 
This allows the part of the supersymmetry algebra that closes on the hamiltonian to 
be realized exactly. A wealth of formulations and studies along these lines can be 
found in the literature. As examples, I cite [40,43,93-100]. See also references in [76]. 
For d > 2 field theories, hamiltonian formulations have traditionally met with only 
limited success. 

3 Deconstructed (2,2) 2d SYM 

In this section one of the deconstruction models of CKKU will be discussed [24]. To 
lay a groundwork for the discussion of the lattice theory, I begin with a brief review 
of the continuum, or, target theory. For the most part this discussion will be classical, 
merely defining the action and the fields that appear in it. 

3.1 The continuum theory 

The theory that I discuss is (2,2) 2d SYM. The "(2,2)" denotes the supersymmetry of 
the theory, generated by supercharges Q±,Q±. The supercharges comprise two left- 
handed 2d Weyl spinors, Q^,Q^, and two right-handed 2d Weyl spinors, 
If equations of motion are imposed, one finds that the notation "(2,2)" equivalently 
denotes the number of left-moving (depending ont + x) and right-moving (depending 
on t — x) 2d supercurrents; see Appendix El The supercharges are derived from the 
supercurrents in the usual way. Often (2,2) supersymmetry is denoted as A/" = 2 
supersymmetry. However, I prefer the former notation since it is possible to build 
chiral models; e.g., models with (0,2) supersymmetry. As an example, generic lOd 
—>■ 4d heterotic sting compactifications with A/" = 1 4d target space supersymme- 
try have (0,2) worldsheet supersymmetry. Compactifications with (2,2) worldsheet 
supersymmetry are the exception, corresponding to so-called standard embeddings. 

The (2,2) 2d SYM theory is most easily obtained as the dimensional reduction 
of A/" = 1 4d SYM to 2d. Note that A/" = 1 4d SYM is nothing but 4d YM with a 
single massless adjoint Majorana fermion. It is minimal, in the sense that the (four) 
supercharges organize themselves into just one 4d Majorana spinor. 4d theories of 
extended supersymmetry, where Af > 1, have more than four supercharges, corre- 
sponding to more than one Majorana spinor. The dimensional reduction to 2d is the 
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naive one: compactify two dimensions on a torus, and truncate to just the zeromodes 
in those dimensions. Equivalently, impose that fields only depend on two coordinates 
t,x. Then 

\ [ dtdx(fy{- • • ) ^ / dtdx{- • • )> (3-1) 

94. J 94: J 

where V2 is the volume of the 2-torus. Consequently the 2d coupling is = g\l'V2- 
It follows that [(^2] = 1, where [■■■] indicates mass dimensions throughout this re- 
port. From the 2d point of view, the two vector-boson modes along the compactified 
dimensions are considered as real scalars in the adjoint representation. Thus, the 4d 
field-strength decomposes as 

Fixv ^ {Fij-, ^i^aA^a.s^-, /i, = 0, 1, 2, 3, i = 0, 1, a, 6 =1,2, (3.2) 

where Sa are the scalars and DiSa = diSa + i[vi,Sa]- Additionally, it is convenient to 
define the complex scalar s = (si -|- is2)/V^- 
The (real-time) action is: 

S = I (fx ^Tr -{Ds) ■ {Ds^ + V2iip_D+ip_ + y/2iip+D_ip+ -\f ■ F 
J 92 . 4 



(3.3) 



Here, ■0± are 2d Weyl spinors, which are 1-component objects. Also, D± — d± -\- 
^b±,-]) ^± — (9t ± dx)/V2 and v± — (vt ± Vx)/V2. The field strength, Fqi — 
dtVx — dxVt + i[vt,Vx\ = E, is just the color-electric field, since we are in 2d. This 
action contains a chiral U{1)r symmetry:''^^ t/jL — > e^'^'^L, fpR — e'^'^t/jR, s e'^'^'^s. 
The (2,2) supersymmetry of the action is 

6v. = tV2{ij+^+ - e+V'+), Sv+ = tV2{^-^_ - ^^-), 
Ss = iV2{^M - = ^V2(V^+e- - e+V'-), 

= 2^_D+s^ - + i[s\ s]), (3.4) 



^°The subscript "R" is traditional, denoting a flavor symmetry that does not commute with 
supersymmetry. 
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where ^±,^± are infinitesmal Grassmann parameters. This is an invariance provided 
the equations of motion are imposed (on-shell supersymmetry) . If auxihary fields are 
introduced, one can avoid the use of the equations of motion (off-shell supersymme- 
try). To obtain the on-shell formulation, one must impose the so-called Wess-Zumino 
gauge. This field redefinition eliminates auxiliary fields but leaves the usual gauge 
symmetry unbroken. The on-shell supersymmetry then involves a mixture the off- 
shell supersymmetry and a gauge transformation. (For further details, the reader is 
referred to [4].) 

After continuation to Euclidean space, a slight change in notation, and the intro- 
duction of 2d Dirac 7-matrices, the action can be written in the form: 



S = I (fx — ^ Tr 

9i 



+iV2{Ms,^pR]+MsKi^L]) + \[s\s]^ . (3.5) 
This imaginary-time formulation is the target of the lattice theory. 



3.2 Construction of the daughter theory 

As summarized briefly in ^2.11 the method of building the class of models described 
here is based on deconstruction of extra dimensions [21,22]. These lattice construc- 
tions are all arrived at by orbifold projections of supersymmetric matrix models; 
i.e., in each case we quotient the matrix model by some discrete symmetry group 
of the theory. Degrees of freedom that are not invariant with respect to the orb- 
ifold generators are projected out. In the case of the Euclidean spacetime lattice 
formulations [24,26,30], the matrix models are Od: they are obtained by reducing 
all dimensions of a Euclidean field theory. In the case of the Hamiltonian (spatial) 
lattice formulations [23], the matrix models are Id: the matrix partition function is 
a functional integral corresponding to a quantum mechanical system with continuous 
time. 

Here I will focus on the orbifolded supersymmetric matrix models that define a 
Euclidean lattice formulation of (2,2) 2d super- Yang-Mills with U{k) gauge group. 
The mother theory is a nonorbifolded U{kN'^) supersymmetric Od matrix model. The 
daughter theory is obtained by orbifolding the mother theory by a x Zn symmetry 
group. This leaves intact (among other things) a U{k)^^ symmetry group. This can 
be associated with a symmetry that acts as an independent U{k) at each site of an 
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N X N lattice. Thus, we obtain an N x N lattice theory with U{k) gauge symmetry. 
Classically, the continuum theory is obtained by studying the daughter theory about 
a particular minimum of the bosonic lattice action. This aspect, and its extension to 
the quantum theory, will be discussed in more detail in ^3.31 below. 
The action of the mother theory is just 

= ^ Tr (^VmnVmn + i'^m[Vm, V^]^ , m = 0, . . . , 3, (3.6) 

where Vmn = i[vm,Vn], and cr^ = are the Euclidean Weyl matrices, with 

CTj {i = 1, 2, 3) Pauli matrices; later we will also make use of the conjugate Weyl 
matrices am = (1, — ^cTj). Each of the quantities appearing in ()3.(jj) is decomposed on 
a basis of hermitain matrices of order kN"^: 

t;^ = CT^ ^ = V'"T", ^ = rT'', a = 0,...,{kN'y-l. (3.7) 

This structure is inherited from the U{kN'^) YM gauge theory, with adjoint Majorana 
fermions, from which the reduction ()3.6|1 springs. It is assumed that we work in a 
basis where Tr T"T^ = 5"^. 

Symmetries play an important role in the construction of the lattice theory. Since 
the mother theory is obtained from a 4d — Od reduction, all symmetries become 
global. Thus, the remnant of gauge symmetry is just 

Vm U^mU, tfj U^^U, tfj U^U, U G U{kN^), (3.8) 

where the matrices U are taken in the kN^ x kN"^ defining representation. 

Since the gauge symmetry is now global, no problems with infinities arise from 
gauge orbits (at finite A^). Thus, in the mother theory partition function, there 
is no need for gauge fixing. The exact value of the SU{kN'^) part of the mother 
theory partition function is known [101, 102]. It is finite — with a value given by 
the D = A,N —>■ kN"^ case of Eqs. (26)-(27) of [101] — in spite of a noncompact 
classical moduli space. (Entropic effects render the quantum moduli space compact; 
cf. App. |B]) The t/(l)diag part of the mother theory partition function is ill-defined, 
as it is a massless free theory. It is completely decoupled from the SU{kN^) part, 
and is irrelevant to all the considerations that follow, though one must be careful to 
factor it out, as will be discussed. 

The model also possesses an SU{2)'^ x U{1) symmetry that is inherited from global 
symmetries of the 4d theory. The symmetry transformation is: 

a^Vm LamVmR\ i' ^pL'^e~'\ i) e'^^Ri), L,Re SU{2). (3.9) 
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This corresponds to the ?7(1)_r and the Euchdean rotation group S'0(4) ~ SU{2)^ 
of the 4d SYM. Finally one has the Od reduction of supersymmetry. This I write in 
terms of infinitesmal 2-component Grassmann parameters k, R: 

= -^K{amOn - On(Jm)Vmn- (3.10) 

Since the mother theory is formulated without auxiliary fields, it is necessary to use 
the equations of motion to obtain the supersymmetric invariance. In this Od case, 
these are just (graded) algebraic constraints: 

= [Vra^Vmn] - l^T {^'a, i^a} = ^m[Vm,i^] = [Vm,i']^m- (3.11) 

In the partition function of the matrix model, ip and ip are independent integration 
variables, not necessarily related by a conjugation constraint. Consequently, the 
supersymmetry transformation parameters k and R that appear in ()3.10p can also 
be taken independently. E.g., there is (formally, at least) nothing inconsistent with 
K 7^ 0, R, = 0. It is still a symmetry of Z. These statements are important below 
when I discuss supersymmetry after orbifold projection. 

Also important below will be the decomposition U{kN'^) C U{k)^^ . I write the 
(hermitian) generators of U{kN'^) in this basis: 

(T°)f,^^^ryy^i', fl,U,fl',u' = 1,. . . ,N, = 1, . . . ,k, 

a = 0,...,{kN'y-l, T° = -^W, TrT"^o = 0. (3.12) 

CKKU exploit a Zj^ subgroup of the SU{2y x U (1) in for the orbifold action, 
and embed it into a subgroup of U{kN'^) symmetry ()3.8j) to break U{kN'^) U{k)^^ . 
The Zli action is generated by exp{2TTira/N) with A^-alities (a = 1,2) that I now 
describe. 

First we define the combinations 

_ vo-tV3 _ -V2-ivi , . 



and fermion components 

7/; = I 



V=(^), i' = ia,P). (3.14) 
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X 


y 


A 




a 




ri 


1 








-1 


1 





r2 





1 





-1 





1 



Table 1: A^-alities of the fields w.r.t. the Z"^ subgroup that CKKU select from 
SU{2f X f/(l). 

From (j3.13p . one finds that 

V = VmO-m = Vo +iv-(7 = V2 \ (3.15) 

\y^ X J 

With these expressions it is straightforward to work out the action ()3.(i|l in terms of 
the new variables. 

The A^-alities of these fields are given in Tabled They are generated by a diagonal 
subalgebra of the su{2)l © su{2)ji © u{1)y associated with ()3.9p : 

n = ~L^ + R,-Y, r2 = L^ + Rs-Y, (3.16) 

as can be checked. Note that the fermion A is neutral. It turns out that in the CKKU 
approach the number of fermion fields that are neutral w.r.t. the orbifold group is 
identical to the number of supersymmetries that are left intact. This simple rubric 
illustrates the usefulness of the orbifold technique. 

The action is embedded in the gauge group as follows. We write 

fi = diag(w,cj^...,tu^), u = exp{2m/N). (3.17) 

Then define 

Ci = fi©l7v©lfc, C2 = Ijv ® fi® Ife- (3.18) 

Thus Ci distinguishes only the indices /x, /i' in ()3.12p and C2 distinguishes only the 
indices z/, u' in ()3.12|) . The action of the orbifold group (with embedding) on any of 
the fields in Table [T] is (a = 1, 2): 

$ ^ e^^'''^/^Ca<l>Cl. (3.19) 

Keeping only the components that are neutral under this orbifold action corresponds 
to the "orbifolded" theory, which will be referred to as the daughter theory, in keeping 
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with the appelations of CKKU [24]. I now examine the exphcit form of the daughter 
theory action. 

It is useful to be very exphcit about ()3.19p . Using the notation of ()3.12p . each field 
in Table HI can be written in the form ^ ^^u^i-^' y ^i' = ^'^{T")^ u^i.yy^ii. Then ()3.19|) is, 
explicitly, 

: $;.,.,i;;.',.',.'^e2-('-^+^-^')/^$^,,,,;^,,,,.. (3.20) 

Thus the fields that are retained after the projection are just those that satisfy 

ri+ fi- fi' = mod A^, ra + z/ - z/' = mod N. (3.21) 

For instance, from Table ^ we see that components of x must satisfy 

XtM,u,i;ii'y,i' 7^ iff jj,' = + l mod A^, u = v. (3.22) 

The constraints ()3.21|) lead to a x periodic lattice interpretation of the 
surviving components of ^ ^,u,iy y ,i' ■ For instance, ()3.22|1 leaves just components 
that can be interpreted as link fields that stretch^^ from the site m = (yU, v) to 
m + i = (/i + 1, z/). Similarly it is easy to convince oneself that y yields link fields 
stretching from m = (/i, z/) to m + j = (/i, z/ + 1). Note also that x'^ and are 
link fields oriented in the opposite direction. By inspection, a yields link fermions in 
the i direction, (3 link fermions in the j direction, and ^ link fermions in the —i — j 
direction. Since fJ^ = fi' and v = v' for A, it is to be interpreted as a site fermion. It is 
convenient to decompose the fields into k x k matrices labeled by site indices of the 
A^ X A^ lattice: 

•^^,u,i;fj.+l,i/,j = {•^(p,i^)^i,j = (•^m)i,j) CtC. (3.23) 

I summarize these conclusions in Table El 

The U{k)^^ subgroup of the gauge symmetry ()3.8|1 that survives is just 

Mj^Xm-Um+i, '^ml/m'Wm+i) G^C. (3.24) 

as dictated by Table |21 Here, the U{k) fundamental representation matrices. 

The fields Xm, T/m, etc., may written in terms of a Hermitian basis of k x k generators 
(/X = 0, . . . , fc^ — 1 and a = 1, . . . , A;^ — 1): 

T^G |y|lfc,T»|, (T'^)t = T^ x^ = x^^,T^ |/^ = C^^ etc. (3.25) 



11 



Note that i = (1, 0) and j = (0, 1) 
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field 




•^m yln. '^rn /^m 




orientation 


* i 


-i -J i 3 


—I — J site 



Table 2: Classification of components that survive, in terms of lattice fields. Fields 
with negative orientation point back to site m, rather than starting from m. 



I define the T°- such that Tr {T^T'^) = 26'^". Furthermore I define 

Tr {T^T^T") = f""", (3.26) 

and remark that (underlining implies all permutations are to be taken): 

^MPO _ ^/.p_ (3 27) 

Note that only A has components that survive for which fi = n' and u = v' . 
Thus it is the only field that has a nonvanishing projection onto the diagonal Cartan 
subalgebra f/(l)'^^^ of U{kN'^): the diagonal components {i = i') of these site fields 
Am correspond to this subalgebra of the mother theory. In particular, the U{l)diag = 
diag ^i^U{l)i^rn projection survives: 

A"T" 3 A°T° = -^A°l,^., A° = ^ y A^. (3.28) 

^ m 

From the commutator form of the fermion action in ()3.6|) . the A° component is an 
ever-present zeromode; i.e., it vanishes identically from the action and the integral 
/ (iA° would cause the partition function to vanish if were included. Since it survives 
the projection, we are guaranteed to obtain a vanishing fermion determinant, just as 
in the mother theory. This U{l)diag sector is completely decoupled, so there is really 
no reason to have included it in the first place. 

To avoid this zeromode, one could just as well have started with the SU{kN'^) 
mother theory. The only price one pays is that the constraint Tr A = needs to be 
imposed. This is just to say Xl/ii/i ^ii,v,i;ij.,u,i = 0, or equivalently in terms of the site 
fermions. 



5^Tr Am = v^5^A^ = 0, (3.29) 
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where the trace is over the U{k) index. This nonlocal constraint (w.r.t. the lattice) 
leads to a nonlocal fermion action when the constraint is imposed by substitution, 
say, of 

^m={N,N) = ~ -^m- (3.30) 

injL{N,N) 

In App. O I discuss an alternative method that factors out the zero eigenvalue. 
That factorization method leaves the fermion action local at the price of having to 
work in a limit process. For numerical applications, I have found that the cost of 
extrapolating to the limit is negligible. 

After all the projections are taken into account, one obtains the following daughter 
theory action: 

•S* = — ^ ] Tr < ^ (^^m-jS^m-i ~ ^m^ln ~^ ?/m-j^m-J ~ 
y m 

+ V^(Q;ni+J^mZ/m ~ Q^mZ/m+i^m) + V^(/3m^m+J^m ~ /^m+i^m^^m) | • (3.31) 

Taking advantage of the decomposition ()3.25|1 . one can organize the fermionic part 
of action into the following matrix form: 

Sf = -^{o.'L. PL) ■ MZr. ■ f5l • (3.32) 



The fermion matrix Mj^^ is given by (sum over u implied in the entries): 



m,n 



m,n'''m 




_f^iup V _|_ ^P-pv V 
m,ni'm+i ' m,n+jyii 


miny m 


m,n— ji'm 


f/iup fPpi' „u 
''m,n'*'m+j ''m,n+i'^n 



(3.33) 



Here I have introduced the compact notation 

C"n = ^m,nt^"^ (3.34) 

where t'^'^P was defined in (j3.26|) . 

The symmetries that survive have been discussed in [24]. The symmetries of the 
lattice theory forbid non- irrelevant operators that would break supersymmetry in the 
continuum limit. Here I only briefly summarize them. On has of course the lattice 
translation symmetry and U{k) lattice gauge symmetry that are apparent from the 
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daughter theory action. The diagonal U{1)^ subgroup of ()3.9|) clearly commutes with 
the orbifold action and survives. There is a nontrivial Z2 point group that reflects 
the lattice about a diagonal link. On the fields it has the action 

«m Pm', Pm ^ Oim' , Am — ^ Am/, ^m ~^m' , (3.35) 

where m = (mi, 7712) and m' = (m2,mi). Lastly, there is an exact remnant of the 
mother theory supersymmetry. Suppressing site indices, 

6x = —V2iari, 6y = —V^i/Sr], 6X = —i{[x\x] + [y\y])rj, 

6^ = 2^[x^ y^]r], 6x^ = 6y^ = 6a = 6(3 = 0. (3.36) 

As usual, 7] is an infinitesmal Grassmann parameter. It is perhaps disturbing that, 
e.g., X and transform differently. One can relabel as x and treat these formally 
as independent complex matrices, as CKKU have done. This is not necessary, how- 
ever; what matters is that an independent variation of the matrix and its hermitian 
conjugate — i.e., a variation that violates this conjugation constraint — is a symmetry 
of the partition function. This suffices to derive associated WT identities for variables 
that do satisfy the conjugation constraint. 



3.3 Classical approach to the continuum theory 

To obtain a discrete approximation to the continuum theory, one expands about what 
I will call the a- configuration: 

Xm = ^7=1, Vm = ^7^1, Vm, (3.37) 
av2 av2 

keeping g2 = ga and L = Na fixed, treating a as small. (It is easy to see that Sq = 

for this configuration.) That is, we associate a with a lattice spacing, even though 

it arises originally from a specific background field configuration. In this case, one 

finds that the classical continuum limit is nothing but (|3.5|) . as has been described in 

detail by CKKU in [24]. 

Here I merely summarize the map to the continuum action ()3.5|) . At each site we 

define: 

1 1 

X = —={a^^ + si + ivi), y = —{a^-^ + S2 + iv2) . (3.38) 
y2 V 2 



^^It has been shown how to exploit polar decomposition to render gauge variables compact, and 
to simplify the symmetry transformations of the fields associated with the continuum limit [103]. 
For brevity, I do not include these details here. 
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The map to continue variables is then given by: 



S(na) = ■^(■5l,n + iS2,n), fm(na) = (f i,n, f2,n), 

^(na) = (An, ^n)^, ^(na) = ?(«„, 



(3.39) 



The Dirac matrices are: 



7l = (^3, 72 = 0-1, 



73 = cr2- 



(3.40) 



The last is the chirality matrix, and projections are defined with 



Pl 



2(1-73) 



^(1 + 73). 



(3.41) 



As usual, on replaces 



(3.42) 



n 



Also, shifted indices that appear in the lattice action are formally expanded in a 
Taylor series expansion of ()3.39|) . E.g., 



The a ^ limit of the lattice action, once these replacements have been made, defines 
the classical continuum limit. 

3.4 Stabilization of the a- configuration 

The classical lattice action contains many zero-action configurations. As described 
above, CKKU expand the classical lattice action about a particular class of zero-action 
configurations that are characterized by a parameter a, the a-configurations. In the 
limit a — > 0, N ^ 00, the classical action tends to (j3.3j) . Thus a is interpreted as a 
lattice spacing. However, it is dynamical as it has to do with a particular configuration 
for the lattice fields. This strategy is based on the ideas of deconstruction [21,22]. 

The validity of the semi-classical expansion about an a-configuration rests on the 
assumption that it gives a good approximation to the behavior of the full lattice 
theory. But all a-configurations are energetically equivalent. Furthermore, there 
exist other zero action configurations that do not fall into the a-configuration class 




(3.43) 



n=0 
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(shown below). For these reasons, CKKU deform the action by adding an a-dependent 
potential that favors the a-configuration. 

CKKU suggest the following deformation of the bosonic action in an effort to 
stabilize the theory near the a-configuration ()3.37|) : 



Here the strength of the deformation is determined by the quantity /i, which has mass 
dimension 1. It is clear that the configuration ()3.37j) minimizes S'sB- The "continuum 
limit" then includes sending a ^ in this potential. 

Unfortunately, the deformation S'sB breaks the exact supersymmetry of the daugh- 
ter theory action (hence the subscript "SB" = Symmetry Breaking). For this reason 
CKKU demand that the strength of S'sB relative to S'o, conveyed by /i^, be scaled to 
zero in the thermodynamic limit {Na —>■ oo): 



Thus we are interested in the effects of the deformation S'sB subject to this scaling. 

In this circumstance, what becomes of all the pleasing symmetry properties of 
the daughter theory? The import of this question rests on the fact that it is the 
symmetry group of the daughter theory action that guarantees the correct continuum 
limit. Without those symmetries, the effective long distance theory will be more 
generic, and not what we intend. 

According to CKKU, the deformation Ssb is rendered harmless by scaling its rela- 
tive strength to zero in the thermodynamic limit. For this to hold true, it is important 
that the fluctuations about the a-configuration be small and weakly coupled. Support 
for the conclusion that the fluctuations are sufficiently small has been provided by 
the arguments offered by CKKU in §6 of [24]. 

In much of what follows I will specialize to the case of U{2). This is merely because 
it is the simplest case and the most efficient to simulate. In this special 
will be unconstrained 2x2 complex matrices. 

3.5 Renormalization 

I will now sketch the simple argument that has been used by CKKU to demonstrate 
that no non-irrelevant operators that break the continuum symmetries can be gener- 




(3.44) 



(3.45) 



/i = l/cL, c = fixed. 



(3.46) 
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ated in perturbation theory. Afterward, I will examine the merits of the argument in 
more detail. 

Consider the form of the perturbative Wilsonian effective action, Sefj = S + AS. 
The correction AS* will be a sum of monomial operators Op ^ that respect the exact 
lattice symmetries. Here, [Op^a] = P, and a labels different p-dimensional operators. 
Quite generally, 

A5 = V ^ [ dPxde Cp.„Cp,„. (3.47) 

Here, / dO is the superspace integration associated with the = subalgebra that 
is preserved on the lattice. Note that [(72] = 1 and [dO] = 1/2. Hence on dimensional 
grounds [Cp^a\ = (7/2) — p. Suppose that g2 and a are the only dimensionful quan- 
tities in the bare lattice action. (I will return to the additional parameters fi and L 
below.) Since, by assumption, the coefficients Cp^a have been obtained from a purely 
perturbative computation, we can express them as a power series: 

Cp,„ = a^-^5^cp,,., ((?2a)'^ (3.48) 

e 

where Cp^a-/ are dimensionless coefficients that are independent of g2a- Thus the 
counterterm corresponding to Op^a at the ith order in perturbation theory has UV 
degree D = —p + | — 2£. Operators that will survive in the continuum limit thus 
satisfy 

P<l-2i, (3.49) 

corresponding to counterterms with UV degree D > 0. On the other hand, in pertur- 
bation theory Op^a will consist of monomials of the elementary fields and spacetime 
derivatives of them. It follows that p > 0. Hence ()3.49j) can only be satisfied for 
i = 0,1. Recalling that i serves as a loop-counting parameter, we see that only at 
1-loop is it possible to generate counterterms with UV degree D > 0. This is nothing 
but the usual result: a theory for which all couplings have positive mass dimensions 
will be super-renormalizable. 

CKKU then show that the only operators that satisfy this requirement are just 
those that already appear in the target theory action and the constant operator, 
equivalent to vacuum energy. The latter can be removed by a suitable normal-ordering 
prescription. Thus one sees that in perturbation theory the target theory action is 
obtained without fine-tuning. 
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Now I include the parameters fi and L in the argument. Since ^ = c/L with c 
fixed, there is only one new dimensionful parameter, which I take to be L. Then 
fl3.48|) generalizes to: 




(3.50) 



Marginal and relevant operators thus satisfy 

7 

p< --2i + m. (3.51) 

It is assumed in the CKKU arguments that one can safely take L — > oo. One must 
believe that there are no m > terms. What is the justification for this? 

In the a — i> limit, the lattice action differs from the continuum action only by 
0{a) terms; for instance, the higher order derivatives that appear in (j3.43|) . Let us 
rescale the boson fields so that they are dimensionless: Vm — > g2Vm, s — > g2S. Taking 
into account the form of the lattice action, it is easy to see that on dimensional 
grounds the extra terms are 0{adm) and 0{ag2). The latter are already accounted 
for in the last factor of (|3.50|) . and do not give rise to the L/a dependence. The 
0{adm) terms only make a significant contribution in the UV domain of integration, 
and should not affect the IR behavior. It follows that no new IR divergences are 
introduced by the irrelevant 0{a) suppressed operators. For these reasons I expect 
that the IR regulator that is used to define the continuum theory will also protect 
the L —>■ limit of expressions such as (j3.50|) : the renormalizations that are sensitive 
to this IR cutoff are just those that already appear in the continuum theory. 

A more heuristic way to put the argument is as follows. The modifications from 
the lattice action all involve irrelevant operators. These are only important in the 
UV, based on renormalization group arguments. It follows that they will not give rise 
to new IR divergences. This is not rigorous and I feel that it deserves a more careful 
examination that is beyond the scope of this work. 

3.6 The fermion determinant 

In this section I will discuss important properties of the fermion determinant. Upon 
examination of ()3.33|) . one suspects that Sp is not hermitian. Thus one would expect 
that the fermion determinant will have a complex phase that depends on the boson 
^•^That this should be true was suggested to me by Joshua EUiott. 
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background. This suspicion was confirmed in [25], based on numerical studies. The 
phase of the determinant is complex and background dependent, as I now review. 

Let det M be the determinant of the fermion matrix with the ever-present zero- 
mode removed. This can be obtained through the projection ()3.29p described above — 
equivalent to having started with an SU{kN'^) mother theory; or, it can be obtained 
from a limiting process on a deformation of M, as described in App. O Alternatively, 
one can just compute the eigenvalues of M and compute the product of those that are 
nonzero. The latter two methods were used in the numerical studies performed in [25]. 
Then the quantity that needs to be studied in relation to a complex phase or sign 
problem is arg det M, the phase of the product of (generically) nonzero eigenvalues. 

For the U{2) case I take 

T'^e{l2,a'^}, a;^ = x°^ + x>", y^ = yl + y^^a\ (3.52) 
Then the fermion matrix is given in ()3.33|1 . with 

^000 _ ^ooo _ 0, = t"^^ = ie"^^. (3.53) 

The lattice theory is obtained by expansion about the a-configuration: 

^m = ^ + ---, < = ^ + ---- (3.54) 
av2 av2 

In the study of det M that was conducted in [25] , I scanned over a Gaussian boson 
distribution where a;^,?/^ had a nonzero mean l/a^/2 = 1. The remainder of the 
bosons were drawn with mean zero. All bosons were taken from distributions with unit 
variance. Though this is not the same distribution as would be generated dynamically 
in a Monte Carlo simulation, it does give a preliminary survey of the behavior of the 
phase of the fermion determinant. 

For a set of 10^ draws on the bosons of this 2x2 lattice, I binned arg det M over 
its range, with bins of size vr/lOO. In Fig. ^ I show the frequency for each bin, as a 
fraction of the total number of draws. The product of the nonzero eigenvalues has 
arbitrary phase. Consequently, there is potentially a very serious problem with can- 
cellations in any Monte Carlo simulation of this lattice system. On general grounds, 
the cancellations will require an exponentially increasing number of samples to get 
accurate expectation values as the system volume is increased. 
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Figure 1: Frequency distribution F{(j)) for = argdetM, for 10^ random (Gaussian) 
draws, binned into intervals of vr/lOO. The distribution of is seen to be nearly 
uniform. These results are for the U{2) lattice theory, with 2x2 lattice. 

3.7 Comparison to U{k) and SU{k) matrix models 

It has already been mentioned that the nonorbifolded supersymmetric U{k) matrix 
models — a mother theory if k = k'N"^ — contain ever-present zeromode fermions com- 
ing from the diagonal U{1), as can be seen from ()3.6|) . Then det M[/(fc) = 0. However, 
it is easy enough to just work with the nonorbifolded SU (k) matrix model, so that 
■0°, ip^ are never in the theory to begin with. Then with appropriate conventions 
det Msu{k) > 0. A proof of this result has been given in an appendix of [104]. 

One might wonder how the generic phases of Fig. Q are obtained in the orbifolded 
matrix models, given the positivity of the fermion determinant in the nonorbifolded 
SU{k) matrix models. That is: If the fermion determinant of the mother theory is 
positive semi-definite, how can the fermion determinant of the daughter theory have 
non-definite complex phase? 

Firstly, one should note that the proof given in [104] shows only that for each 
eigenvalue c of Msu{k) there also exists an eigenvalue e*, and that these always come 
in pairs. The proof relies essentially on the relation 



(3.55) 
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On the other hand, using ()3.15|) . and ()3.19p as it apphes to x, ?/, it can be shown that 
the Zn X Zn orbifold action on the mother theory bosons is generated by 

a(rpa\ -i-Kci^ IN a -i-Kcfi IN q iT'^\ 

'J V jipujp'u' c c i.Lp^y± Ji/iujfi'u'^^^ip'i 

a(rpa\ ^ i-Ka^lN a -i-Ka^lN q (n-^ci\ 0~1 (1. ^P.^ 

where Q was defined in fl3.17|l . Thus because of the exp{±ina^ /N) factors, the orbifold 
projection does not commute with the operations of the proof given in [104]; i.e., 
Eq. ()3.55p . The fermion matrix of the orbifolded matrix model lacks many of the 
eigenvalues of the mother theory; it should come as no surprise that not all eigenvalues 
are removed in pairs (e, e*). After all, we already know that only some of the zero 
eigenvalues are removed from the U{kN'^) mother theory. For this reason it is not at 
all contradictory that the projected theory has a product of nonzero eigenvalues that 
is not positive, nor real. 

3.8 Source of the complex phase problem 

1 now examine the reasons why the fermion determinant is not real. In this discussion, 
it is important to note that real-time hermiticity takes on a generalized sense in 
the imaginary-time, Euclidean theory. What is needed is the unitary involution G, 
satisfying = 1, that is the continuation of the real-time hermiticity. Background 
for this discussion is provided in App. |Dl 

First consider the Euclidean continuum theory ()3.5|1 . Because of the continu- 
ation to imaginary-time (Euclidean) from real-time (Minkowski), Ie = —itui the 
hermiticity of the real-time action is replaced by a conjugation that takes t —t 
in the imaginary-time formulation. Furthermore, the form of the Dirac adjoint, 
ip = ip^'jo = —'ii''^l2, requires that the conjugation of ip,ip in Euclidean space in- 
volve factors of ^72. Finally, the identification V2 = Wq requires that the conjugation 
change of the sign of V2, since Vq was taken to be hermitian in real-time. 

It is easily checked that the operation of hermitian conjugation for the real-time 
action continues to the following generalized conjugation in Euclidean space: 

6'?/'(x, t) = i'ip{x, -t)72, Qi>{x, t) = i'^2i>{x, -t), 

9fi(x, t) = fi(x, — t), 6f2(x, t) = —V2{x, —t), 

Qs{x, t) = s^{x, -t), Qs\x, t) = s(x, -t). (3.57) 

The operator also conjugates complex numbers and transposes matrices and ex- 
pressions involving fermions. E.g., B(^/'M^) = QijjM^Qijj for a constant matrix M. 
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It is easily checked that 6^ = 1. It is a simple exercise to verify that the Euclidean 
action ()3.5|1 is self-conjugate under 6. That is, 65* = S. This must, of course, be true 
since the real-time action is hermitian and B is just a straightforward continuation 
of this to imaginary-time. 

Now I translate this conjugation to the lattice. I will choose to associate the j 
direction with imaginary-time X2 = t, and the i direction with space Xi = x. Then 
comparing (j3.57|) to (j3.38j) - (j3.39|) . one sees that the simplest map of the continuum 
conjugation^^ is just: 

= Vln'^ ©-^m = -/5m', ©^m = -«m', 
©ttm = — ^m', ©/^m = —Am', 

with m' = (mi, —1112) where m = (mi, ^12). (3.58) 

That y does not transform into y"^ just follows from the f 2 — > — f 2 property of ()3.57p 
and the complex conjugation property of G when acting on complex numbers. 

It is straightforward to check that the lattice action St = S + Ssb, given by ()3.3H) 
and ()3.45|1 . is not self-conjugate with respect to ()3.58j) . On the other hand, as has 
already been mentioned, we can express the approach to the continuum limit in terms 
of a power series of continuum expressions: 

5-^ = + + a25(2) + . . . . (3.59) 

Here, S^^^ is just the continuum action ()3.5|) . As has been stated, e5(°) = It 
follows that 

OSt -St = Q{aS^^^ + a^S^^^ + ■■■) = 0{a). (3.60) 

Thus, the violation of the continuum self-conjugacy is entirely due to the irrelevant, 
0{a) suppressed operators. 

In particular, it is not difficult to show that 

QSf - Sf = X] '^^ I ^ "mVs xj^Am - Pin^j>^in'V 2 yln 

+am+i^mV^ym + /^mV^Xm^mj- (3.61) 



""^^We could imagine generalizations that smear nearby sites in a unitary fashion. 
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Here, forward and backward difference operators are used: 

ciW 2 -^m ~ 3^m+j ^mi '^^2 -^ra ~ ^m—j- (3.62) 

From ()3.61|) one sees that terms that violate self-conjugacy of the fermion action are 
0{pa) when we pass to momentum space. These will not be small for the modes 
with p = 0{a~^). The fermion determinant includes the product of eigenvalues 
associated with these modes, independent of how small we take the lattice spacing a. 
It is therefore not surprising that we do not recover a real fermion determinant. We 
always have eigenvalues that reflect the lack of self-conjugacy of the fermion action. 

Lastly, I remark that the boson lattice action ()3.44|1 is not self-conjugate w.r.t. to 
()3.58|) either. The violation is, according to the general arguments made above, 0{a). 
On the other hand, the boson lattice action is self-conjugate under ordinary hermitian 
conjugation. On the a-configuration, the fermion lattice action reduces to 

= \Y.^'' {(Vr«„. + V2Pm)Xrn + {V+a^ - VtPm)U} ■ (3.63) 

^ m 

Then if we define 

'^m -^m? Cm Q^m Oimj Pmi (3.64) 

one finds that Sp^ = Spo- This property is easily extended to the non-irrelevant 
interaction terms involving the vector-boson and scalars. On the other hand, it can 
be shown that the irrrelevant, 0{a) suppressed terms do not enjoy this property. This 
too is consistent with the fact that the fermion determinant is complex. 



4 Deconstructed lattice (4,4) 2d SYM 

As a further example, I now consider another 2d SYM theory. In [26], CKKU proposed 
a lattice action for (4,4) 2d SYM. In [105], some aspects of the fermion determinant 
were examined, with results very similar to those for the (2,2) model. I will briefly 
present those results below. I will more extensively review my efforts [27] to study a 
rather fundamental aspect of the construction: the proposed emergence of the type 
of dynamical lattice spacing that was described for the (2,2) model above. The semi- 
classical approach lends much support to the assertion that the continuum approach 
is exactly as CKKU have proposed. It would be nice to obtain this same result from 
Monte Carlo similulations of the model, which I have attempted. Unfortunately, the 
complex phase problem has left me with only supportive but inconclusive results, as 
will be seen. 
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4.1 The continuum theory 

One begins with 6d YM minimally coupled to single massless symplectic Majorana 
fermion in the adjoint representation. This is A/" = 1 6d SYM; then (4,4) 2d SYM 
is obtained as its dimensional reduction to 2d. Similar to the reduction described 
above, (/I = g\lR^, and again [(72] = 1. Note that the 4 vector boson modes along 
the compactified dimensions are considered as real scalars from the 2d point of view. 
As usual, the notation "(4,4)" denotes that in the 2d theory there are 4 left-moving, 
left-chirality super currents, and 4 right-moving, right-chirality super cur rents. 

Continuing to a Euclidean formulation, one finds that the 2d Euclidean action is 



S = I (fx — Tr 

92 



+^i[{so5ij + Z73S ■ (Tij),tjjj] - ^[s^, Su]^ 



(4.1) 



where [fi = 0,1,2,3) are hermitian scalars, F is the 2d YM field strength, and 
ipi {i = 1, 2) are 2d Majorana fermions, all in the adjoint representation. For our 
purposes, the gauge group will be U{k). 



4.2 The mother and daughter theories 

The mother theory is obtained from the Euclidean A/" = 1 6d U{kN'^) super- Yang- 
Mills. The action is dimensionally reduced to Od to obtain a [/{kN"^) matrix model. 
The matrix model naturally possesses 5*0(6) ^ SU{A) Euclidean invariance. Next 
we note SU{4:) D U{1) x f/(l) D Zn x Z^. This is very much like SO{A) ~ SU{2) x 
SU{2) D U{1) X U{1) D Zn X Zn of the (2,2) model above. CKKU identify a 
homomorphic embedding of Zn x Zn into the U{kN'^) gauge symmetry group, just 
as in Eqs. ()3.17j) - ()3.18j) of the (2,2) theory. Using this, a Z^ x Zn orbifold projection 
is performed to obtain a U{k)^^ Od quiver (or moose) theory. Due to the similarity 
to the (2,2) model, I will not review the details of this procedure, but just state the 
end result. 
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The bosonic part of the daughter theory action is 
'1 



+2(l/n2n+i - 2;nl/n)(4+jl/n " 1/^4) 

+2(2nXn - Xn2n+i)(x|,4 - 4+i4) (4-2) 

Here, Xm, are bosonic lattice fields that are k x k unconstrained complex 

matrices. The U{k)^^ symmetry is nothing but the local U{k) symmetry of the 
lattice action 5*0, with link bosons x-^i in the i direction, link bosons y-ai in the j 
direction, and sites bosons Zm-, all transforming in the usual manner: 

Canonical mass dimension 1 is assigned to Xm,ym, ^m, whereas g has mass dimen- 
sion 2. 

The fermionic lattice action can be written in the form 



(4.4) 



The elements of the fermion matrix are: 



Wn)l,l 



mn/ J^i^ 



(M, 

(MZ)l,4 
(MZ)2,3 
(MZ)3,1 
(Mrj3,4 



4,2 



(MZ)2,2 



(MZ)3,3 



m,n~n+2 ' ''m,n^n' V^"mn/1,3 



mnM,4 — 0, 
''m,nyn+ 



m,n— r m' 



(MZ)2,1 



(MZ)2,4 



m,nyn 



''m,ni/n ''m.n+ji'in l^^^i 



11111)3,2 



''m,n*''n 



_ fPP" 

m,n+j"n' 

rppu V 

''m,n^n' 

_ +PP'' 

mjii—jyim 

m,n+i'''rn 



+pvp —V I y_p^- . 

''m,n'^n+J ^ ''m,n-i*^m' 



PP \ 
mnj4,3 



''m,n^n+?+J ''m,n^n- 



(4.5) 



The notation is as in the (2,2) model above. Note in particular the definitions ()3.26p 
and dSSl- 
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This matrix has two ever-present fermion zeromodes. The first is associated with 
left multiphcation (the transpose "T" merely indicates a row vector, for consistency 
with (Q): 

(CeCn,eCn,A^;r = (0,0,0,AO^ V n. (4.6) 
The second is associated with right multiplication: 

Km , ' ^3,m , Xii.) = (0 , , ijS'^ , 0) V m (4.7) 

The eigenvalues of M are defined through left multiplication. Only ()4.6|) corresponds 
to an ever-present eigenvalue of M, whereas ()4.7|) corresponds to an ever-present 
eigenvalue of M"^. Because the matrix M is not hermitian, it is diagonalized by 
M ^ D = UMV with U and V independent unitary matrices. When this is done, 
the diagonal matrix D has just 1 zero on the diagonal, corresponding to a column of 
V and a row of f/, from ()4.6|) and ()4.7|) . These ever-present zeromodes originate from 
site fermions that are associated with the t/(l)diag subgroup of the mother theory. 
Having identified these modes in the daughter theory, it is straightforward to project 
them out as in ()3.29j) . 



4.3 Construction of the lattice theory 

To obtain the lattice theory, one proceeds in analogy to the (2,2) theory. An expansion 
is made about the a- configuration 



1 



1 



1^/2 



Vr. 



0, Vm, 



(4.^ 



keeping g2 = ga and L = Na fixed, treating a as small. It is easy to see that Sq = 
for this configuration. As before, we associate a with a lattice spacing, even though it 
arises from a background field configuration. One finds that the classical continuum 
limit is ()4.H) . as has been shown in detail in [26]. 

In order to make the configuration ()4.8|) energetically preferred, CKKU utilize a 
deformation that is just a slight modification of the one for the (2,2) theory, Eq. ()3.45|) : 



Sb 

SsB 



So + Si 



SB, 



2g^ 



(4.9) 
(4.10) 
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That is, a term is added to stabilize the 2;-field about zero. As in the (2,2) case, Ssb 
breaks the exact supersymmetry of the daughter theory action. Once again CKKU 
demand that the strength of Ssb relative to 5*0, conveyed by /i^, be scaled to zero in 
the thermodynamic limit ()3.46p . 

As in the (2,2) case, I will generally specialize to U{2) gauge theory in what 
follows. 

4.4 The semi-classical analysis 

As will be explained below, the bosonic action of the undeformed daughter theory 
satisfies 5*0 > and a vast number of nontrivial solutions to 5*0 = exist, not all of 
which are gauge equivalent. In fact, the space of minimum action configurations, or 
moduli space, is a multi-dimensional noncompact manifold with various branches — 
classes of configurations. In ^4.4. H some of these minima will be illustrated. In ^4.4.21 
I will consider the modifications induced by the deformation Ssb, which has the effect 
of lifting all flat directions in moduli space, modulo gauge equivalences. That is, the 
a-configuration is the absolute minimum of Sb = Sq + Ssb- However, because the 
deformation is scaled to zero in the thermodynamic limit, it is important to keep in 
mind the undeformed moduli space. Through understanding this classical picture, 
naive expections of what will occur in the quantum theory, based on energetics, 
can be formulated. In §4.4.4| the effect of the fermion determinant det M will be 
considered. It will be shown that the minima of the undeformed bosonic action 5*0 also 
have detM = 0. Significantly, this is also true of the a-configuration. Fluctuations 
about these minima of 5*0 cause the determinant to be nonvanishing; however, for 
small fiuctuations I find that the magnitude of the determinant differs significantly 
depending on which minimum they perturb. This gives some hints as to the relative 
contribution of these various configurations. My conclusion will be that the continuum 
limit is the target theory, but that a renormalization of the dynamical lattice spacing 
occurs. 

4.4.1 Undeformed theory 

Here I neglect Ssb and examine the minima of 5*0. Note that ()4.2p is a sum of terms of 
the form Tr AA'' (the first line involves squares of hermitian matrices). Thus 5*0 > 
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with 5*0 = iff the following equations hold true: 

•^nUn+i Z/n-^n+J Vn^n+j ^nZ/n ^n'^n -^n^n+i 

The set of solutions is the classical moduli space of the undeformed theory. 

Zeromode branch. To begin a study of the moduli space, I isolate the zero 
momentum modes: = x Vn, etc. Then Eqs. 1)4.111) and ()4.12|) reduce to 

[x\x] + [y\y] + [z\z] = 0, 

[x,y] = [y,z] = {z,x]=0. (4.13) 

Eqs. ()4.13p may be recognized as nothing but the D-flatness and F-flatness constraints 
that describe the moduli space associated with the classical vacuum of A/" = 4 4d 
super- Yang-Mills. The equations are invariant with respect to the global gauge 
transformation 

X axa\ y aya\ z — ^ aza^ . (4-14) 

Then it is well-known that solutions to ()4.13j) consist of x, z that lie in a Cartan 
subalgebra of U (k); the proof is reviewed in App.|E| The global gauge transformations 
fl4.14|) allow one to change to a basis where this Cartan subalgebra has a diagonal 
realization. Thus one can think of the moduli space as the set of all possible diagonal 
matrices x, y, z, and all global gauge transformations ()4.14|1 of this set.^^ 

In particular, the zeromode moduli space of the undeformed U{2) theory is com- 
pletely described by 

x = x° + xV^, y = y^ + y^a^, z = z° + z^a^, (4.15) 

with arbitrary complex numbers x^,x^,y^,y^,z^,z^, together with U{2) transforma- 
tions of these solutions. 

Eqs. ()4.11)1 and ()4.12j) also have non-zeromode solutions. I do not attempt to 
present an exhaustive account of them. I will merely point out a few such branches 

-"^^I thank Erich Poppitz for pointmg this out to me, as weU as the branch of moduh space (|4.16|l 
given befow. 

Actually, the gauge transformations 1)4. 3|l that follow from (|4.11|) - (|4.12|) also leave the equations 
(|4.13() invariant. This generates an inhomogeneous solution to 5*0 — 0, and will be important in 
2331 below. 



40 



(4.11) 
(4.12) 



in order to illustrate that the undeformed theory has a very complicated and large 
set of 5*0 = configurations. 

Xm = Z/m = branch. This is the branch of moduli space described by 

Xm = ym = 0, z^ = zl^ + zj^a^, Vm. (4.16) 

Again, z'^, are arbitrary complex numbers. Furthermore, z^ is a site variable and 
thus transforms independently at each site as 

Zm am^^mttL- (4-17) 

It can be seen that this branch affords a vast number of solutions to ()4.1H) and 
(|4.12p : there are iV^ such solutions, modulo choices for z'^, & C and gauge equiv- 
alences. 

Zn = branch. Another branch in moduli space is the following. First I set 
Zn = Vn, and introduce Fourier space variables 

= ^ Vn = ^ JZ^'^ V, uj = exp{2m/N), (4.18) 

k k 

where k = {ki, and fci, k2 G [0, 1, . . . , — 1]. Then taking into account = 0, 
the conditions 1)4.111) and ()4.12|) are equivalent to: 

k 

= J2i'''''^'^''^f^9-^.-i-^'''9-^.^ih), (4.19) 

k 

for all £ = (£i, £2) and £1, £2 G [0, 1, . . . , — 1]. Next I turn off all modes except one 
for both /k and g\s_: 

/k = f^k.k'/k', 9k = 4 -k'fi'-k'- (4.20) 

Here and below, no sum over k' is implied. When substituted into (|4.19|) . only 2 
nontrivial conditions survive: 

= [fl, /k'] + [9\„9-i.'], = /k'(7-k' - uj^'^'^-'^'g^^^h,. (4.21) 

For the U{2) case, solutions exist if uj^'^+S)-^' = -|-]^_ 

For have solutions for /k',fi'-k' diagonal matrices. This includes 

the zeromode solution (k' = 0) ()4.15j) with z = 0. The a-configuration ()4.8|1 is a 
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special case. I want to emphasize that any statement that holds generally for the 
Zn = branch will hold in particular for the a- configuration. 

In the case of lu^'^'^^^'^' = —1 it is easy to see that there are solutions of the form 

/k' = ^_k' = Zg{a^ + ba^), Zf, Zg e C, b e R. (4.22) 

There are many values of k' for which u^'^^^^'^' = ±1. For N even, these are 

3N 

k[ + K = 0,-,N,—. (4.23) 
For odd, k[ + k2 = 0,N are allowed and in the cases where 

,4.4) 

()4.22|1 yield approximate solutions to ()4.2H1 . with an error of order 

Thus in the N oo limit the number of 5*0 = configurations in the z^ = 
class is vast; in fact, it is easy to check that the number of such configurations is 
approximately 2N, modulo gauge equivalences and various choices for the constants 
such as in (KT^ . 

4.4.2 Deformed theory 

Now I consider the supersymmetry breaking deformation S'sB introduced by CKKU. 
To see its effect it is handy to rewrite the quantities that appear in it. Recall that 
complex 2x2 matrix. Dropping the subscript, we can always define 



X 



Then it is straightforward to work out (yU = 0, . . . , 3) 



XX 



a = 



(4.25) 



(4.26) 



Note that 0^'^ are real, and that 0^'*^ is positive definite. With similar definitions for 
(/)^, (/)^, the CKKU deformation is 
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2^ 
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m 
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2^ 
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,2 



(4.27) 
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It can be seen that the deformation drives (j)^, 4>^, (f)^ toward the origin, and 
0m°5 (f'^ toward l/2a^. When (j)^ = 0, it is easy to see that (f)^ = identically. 

To continue the analysis, it is convenient to rescale to dimensionless quantities — 
denoted by a "hat" — using the parameter a: 

g = ga^, ft = fia, 0^ = a^cpl^, Xm = axm, etc. (4.28) 

Then Ssb = 



A! 

9' 



m 



(4.29) 



For any value of the lattice spacing a, the minimum of ^sb is obtained iff 



4>'J = ^'m' = l. 0^° = 0r = 0r = O, Vm. (4.30) 
For example, from ()4.26|) one sees that the conditions involving are just 

Let us examine what additional constraint this places on classical solutions to = 0, 
beyond the restrictions ()4.1H1 - ()4.12|1 of the undeformed theory. 

First 1 note that most of the non-zeromode configurations discussed in §4.4.11 
above are not minima of ^sb. Non-zeromode cases that are minima will be discussed 
in ^4.4.3l below. For the zeromode configurations ()4.15|) . Eqs. ()4.3ip imply that ()4.15|) 
is restricted to the form 

x = ^ diag (e*'^^ e"*'^^) (4.32) 
V2 

and symmetry transformations of this. That is, x is restricted to be an element of the 
maximal abelian subgroup U{lY of U{2), up to an overall factor of 1/V2. Similarly, 
we have for y, 

diag (e*'^^e-*'^'') . (4.33) 



^/2 



Finally, 0^ = implies z'^^ = 0, which has the unique solution = 0, Vm. 

There is a U{1)^ global symmetry of the daughter theory, described in [26], that 
can be used to remove the phases 'jx,ly Apart from global obstructions that are 
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essentially Polyakov loops in the i or j directions, the phases (fx, in the configu- 
ration ()4.H2j) and ()4.HH|1 can be gauged away. It is straightforward to verify that the 
required gauge transformation is ()4.3p with 

am^^ra2 = diag (e^(™i^-+"^2^«)^ g-iCmiv'.+^a^.)) _ ^4_34) 

This sets all Xm,ym to 1/V2 except at the boundaries — the location of which is a 
matter of convention due to translation invariance and periodic boundary conditions: 

XN,m, = -^diag (e*^'^^e-^^'^^), Vma; 

Vn^uN = ^diag (e*^'^^e-^^'^-), Wm,. (4.35) 

For most purposes, I do not expect such vacua to distinguish themselves from the 
trivial vacua in the thermodynamic limit. In any case, global features such as these are 
typical of classical vacua of other lattice Yang-Mills formulations, such as the Wilson 
action. In my simulation study (cf. ^4.5p . I have avoided this issue by restricting my 
attention to the expectation value of a quantity that is independent of these angles. 

4.4.3 Inhomogeneous minima of the deformed theory 

Consider the = branch with uj^^^^^'^ = 1, where I have dropped the prime on k'. 
Following the details given above, one has 

Xm = ^^^-'"/k, 2/m = (4.36) 

with /k, (7_k diagonal matrices. That is, 

/k = 6 + da^ g_^ = b' + d'a\ b,d,b' ,d' e C. (4.37) 
In the notation of ()4.26|) and ()4.28|) . it is easy to check that 

= ^M + MH' '^^''^ = ^,ibd + bd)S^\ (4.38) 

and similar expressions for 0^. From ()4.29|1 one sees that a minimum of Ssb is 
obtained provided 

bd + bd = b'd' + b'd' = 0, + |rf|2 = |5f + |rf'|2 = iL^. (4.39) 
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A simple parameterization of these minima is the following: 

6 = ^e*^cosr/, d = i^e''^smri b + da^ = ^e''^ e'"''', (4.40) 

with similar definitions for the primed variables. Equivalently, 

= ^cu'^-'^e'^e^'"^', ym = ^uj-''-"'e'''e''''''\ (4.41) 
av2 av2 

The resulting 2d x 2d space of minima of the total bosonic action, Sb = So + Ssb, 
occurs for each value of k s.t. 

k2 = -ki mod N. (4.42) 

Next I note that with ()4.39p imposed, 

2 _ _ 

Tr a^x^xl, = —{bd + bd) = 0. (4.43) 

This raises the possibity that ()4.4H) are gauge-equivalent to the a-configuration, or 
more generally the configurations ()4.32|) - ()4.33|) that differ from the a-configuration 
only by edge effects. I now show that this is the case. 

Since by assumption the matrices appearing in ()4.36|) are diagonal, we can restrict 
to the diagonal f/(l)^ subgroup of the U{2) gauge transformations ()4.3|1 . Furthermore, 
the m dependence in ()4.3fij) is an overall phase, so we can restrict further to the 
f^(l)diag subgroup. Thus I take = exp{i(^) in ()4.3|) . The assumption that ()4.4ip 
are gauge transformations of ()4.32p - ()4.33p implies: 

6' = 7x, V = ^x, 0' = 7j^, T]' = ipy, (AAA) 

g«(Cm-Cm+f) _ i^km^ g«(Cm-Cm+j) _ QJ-^"^ (4.45) 

It is not hard to solve the recursion relations ()4.45p : 

C(mi,m2) = C(i,i) - ^ [("^1 - + {2 -mi- ms)] , (4.46) 

with the phase C(i,i) site (1, 1) arbitrary. To prove that these satisfy ()4.45|) . it is 
important to keep ()4.42|1 in mind. 

Similarly, Ssb = solutions exist on the = — 1 part of the = branch. 

From (jOS), 

Xmxl = ^k/l', ymvl = + b'). (4.47) 
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Choosing \zf\'^ = \zg\'^{l + b"^) = N'^ /2o? gives the minima. However, these too are 
gauge equivalent to the configurations ()4.H2j) - ()4.HHj] . The gauge transformations are 
more comphcated due to the nondiagonal nature of ()4.22p and the constraint 

N 

k2 = -ki + — mod A^. (4.48) 

The conditions that must be satisfied by the gauge transformation parameters ()4.3p 
are: 

ttm+i = ^""^ ™ttm (^Cr^), Om+j = ^'""ttm [^(cOS</)Cr"^ + sin(/)Cr^)] (4.49) 

where = tan^^ h with h defined in ()4.22j) . These recursion relations are solved by 
am = e^'^'"a(i,i)[i(cos0a^ + sin0a2)]"^-i(z(T=^)'"i-\ 

TT 

Cm = C(i,i) + ^["^2("^2 - 2mi + 1)A;2 

+ (2(m2 - 1) -mi(mi - l))A;i]. (4.50) 

To summarize, I find that with the addition of the deformation 5*5^, the a- 
configuration is the only one that persists as a classical minimum, modulo gauge 
equivalent configurations and boundary effects. The = branch consists of the 
zeromode branch ()4.15j) with ^ = 0, and gauge transformations of this. 

However, the strength /i of the deformation is scaled to zero in the thermodynamic 
limit. One might wonder if the Xm = Vm = minima of Sq that were lifted by Ssb 
become important again in this limit. I next show that this is not the case, due to a 
suppression from the fermion determinant.^^ 

4.4.4 Determinant on moduli space 

I now show that for each of the boson configurations discussed in §4.4.11 above — 
branches of the undeformed theory moduli space — additional fermion zeromodes ap- 
pear. It will be seen that this is also true for the a-configuration. In configuration 
space, the weights of the immediate neighborhoods of these 5*0 = points are damped 
by fermion effects, since detM ~ 0. (Recall that in M the ever-present fermion ze- 
romode has been projected out.) Of course, even in the continuum theory the trivial 
vacuum — equivalent to the a-configuration — has a vanishing fermion determinant (for 
periodic boundary conditions). What is important is the relative weight of the various 
-"^^This is to be contrasted with what happens if the fermions are not included [106]. 
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neighborhoods of the different minima of Sq- Since the strength of Ssb is scaled to 
zero in the thermodynamic hmit, considerations other than the exp( — Ssb) suppres- 
sion come into play when one attempts to determine the correct configurations that 
should be used for a semiclassical expansion. The determination of the dominant sad- 
dlepoint of the effective action — about which the continuum limit is best defined — is 
a subtle question; it involves an interplay between the size of det M x exp(— 5*^) and 
entropic effects. It will be seen that the a-configuration dominates. 

The a-configuration. It is a simple matter to work out the fermion matrix for 
the a-configuration ()4.8|1 from ()4.5j) . Then one can check that in ()4.4j) 

(ef,n, e£„, ef,n, Kf = 6, 6, A)^ Vn, (4.51) 

is a zeromode for any values of ^j, A. It follows from the considerations of App. [0 
that detMe(f) = det(M(f) + e) oc e^, where v is the a-configuration. Moreover, if 
1 ^ 5 ^ e, we have det M^{v + 0{6)) oc e6^, since the 0(e) eigenvalues not associated 
with ever-present zeromodes are shifted by 0{6). The conclusion of this is that for 
order 6 fluctuations about the a-configuration, the suppression due to approximate 
zeromodes is det M = 0{6^). Next I will compare this to the other points of the 
^0 = moduli space that were discussed in ^4.4. 1[ 

= Vm = branch. Here it is easily checked that M acting on the vector 

(ef,n, eCn, eCn, Kf = (0, 0, 0, XJ^Y (4.52) 

vanishes for any choice of the A^^ Grassmann variables An- This gives A^^ fermion 
zeromodes associated with left multiplication. Thus det M^{v) oc e^^, where v denotes 
the selected point along the Xm = Vm = branch. Again introducing fluctuations 
of order 6, with 1 ^ 5 ^ e, we see that detM^ oc e5^^~"\ leading to a det M oc 
^N^-i suppression. Thus the inclusion of fermions dramatically changes the weight 
of configuration space in the neighborhood of Xm = Vm = branch in the integration 
measure. Relative to the a-configuration, there is an 0{N'^6^^~'^) suppression, where 
the factor A^^ comes from the number of solutions ()4.1(i|l . In the thermodynamic limit 
(A^ oo), this branch has vanishing weight in comparison to the a-configuration. In 
a semi-classical analysis of the continuum limit it need not be included. 

-2m = branch. Due to the gauge equivalence that was established in ^4.4.3| 
we already know that the number of fermion zeromodes for Ssb = configurations 
will be the same as for the a-configuration. However, the configurations with Ssb 7^ 
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also have the same number of fermion zeromodes. This is not surprising, since other 
configurations on the = branch correspond to minima of Ssb with some other 
value for the parameter a. 

First consider the case where _ 2. Then set 

(et„, e£„, eCn, Kf = ^-'^■"6, 6, A)^ V n, (4.53) 

with C,i, A arbitrary Grassmann numbers. It is easy to check that these are fermion 
zeromodes on this branch. The A direction is just the ever-present zeromode. The 3 
directions are particular to this branch. It can be seen that the a-configuration zero- 
modes ()4.5H) are a special case of ()4.53|) . corresponding to k = 0. Clearly, the supres- 
sion of fiuctuations about this branch, due to approximate fermion zermodes, is of the 
same order for all points on the branch: det oc e^^, and hence det M oc 6^. It is also 
simple to see that the fermion zeromodes are related to those of the a-configuration 
by the gauge transformation ()4.45|l - ()4.4(i|l that relates the two configurations, as must 
be the case. 

Now consider the case where cu'^i+'^a — _i Then set 

(etn,C,eCn,A^;)^ = 5^°(a;-'^"6,^'"6,(-)"^-^"^6,A)^ V n. (4.54) 

with ^i, A arbitrary Grassmann numbers. It is not hard to check that these are fermion 
zeromodes on this branch. 

Summary. Due to the gauge equivalence of the a-configuration and other 
So = Ssb = configurations, and due to the suppression of the Xm = Z/m = branch, 
semi-classical expansion about the a-configuration suffices to obtain reliable estimates 
in the lattice theory. In particular, one is interested in the effective lattice spacing, 

4- = (Tr a;ma:L) = (Tr ?/m2/L), (4-55) 

in the a ^ 0, — *^ cxo limit. A renormalization of the lattice spacing has already 
been studied perturbatively about the a-configuration in [107]. Another approach to 
studying this question is by Monte Carlo simulations, which I now take up. 

4.5 Simulation analysis of the dynamic lattice spacing 

In the present subsection I review Monte Carlo simulations [27] that investigated the 
stability of the a-configuration. Some idea of the effect of the fermion determinant 
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has already been gained from the semi-classical analysis of ^4.4.41 The intent of the 
simulation analysis was to pursue a more systematic and complete analysis. I have 
studied 

(0Jn°) = (S^'JL) = Tr (wL)^ = alJ2a' (4.56) 

in my simulations. This expectation value is to be compared to the classical predic- 
tion doni). 

4.5.1 Scaling 

I study ()4.56p along a naive scaling trajectory: 

g2 = a~^g{a) = fixed. (4-57) 

That is, I hold the bare coupling in physical units, g2, fixed. The dimensionless bare 
coupling g is then a function of a that vanishes linearly with a as the UV cutoff is 
removed. 

With regard to /i I follow the instructions of CKKU: I send the dimensionless 
coefficient fi of the deformation S^b to zero as while increasing N: 

fi-^ = cN, c = (4.58) 

where c is a constant. This is equivalent to ()H.4(jj) . 

In the rescaled variables ()4.28p . the coefficient of the undeformed action is l/g"^, 
whereas the coefficient of the deformation is fi'^/g'^, as can be seen from ()4.29|) . The 
relative strength of the deformation vanishes in the thermodynamic limit, when ()4.58|) 
is imposed. 

I perform these scalings for a sequence of decreasing values of a. I then extrapolate 
toward a = to obtain the continuum limit. The length scales (in lattice units) 
are set by g2^ and the system size A^. To keep discretization and finite-size effects 
to a minimum, I would like to take 1 ^ g^^ ^ N, but often violate the bounds 
1 < g^^ < N for specific points where measurements are taken. The reason for this 
is that data outside the optimal window is informative to the extrapolation. 

4.5.2 Problems with the complex phase 

It is essential to include the fermions if we are to draw conclusions for the supersym- 
metric system. The present system suffers from a complex fermion determinant [105]. 
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Figure 2: Average frequency distribution for a = arg det M in the phase-quenched 
distribution for the 2x2 lattice. Error bars are determined by variance in bin counts 
for several blocks of data. 

The integration measure det Me""^-^ is not positive semi-definite and does not provide 
a satisfactory probability measure for Monte Carlo simulations. If we factor out the 
phase e*" = detM/| detM| and use | detM|e~'^^ instead, then we will generate the 
phase- quenched (p.q.) ensemble of boson configurations. Expectation values of an 
operator O in the full theory are formally related to those in the phase-quenched 
theory by the reweighting identity 

(O) = {e'-0)^J{en,.,. (4.59) 

In some lattice systems, the distribution of a in the phase-quenched ensemble is 
sharply peaked. The phase-quenched quantities in the ratio exist and can be measured 
reliably. For example, this is the case in the 4d U{1)l x symmetric Yukawa 

model of [108]. By contrast, the results for the CKKU system are not encouraging. 
The distribution of a = arg det M in the phase-quenched ensemble was found in [105] 
to be essentially flat, as can be seen from Fig. |21 This leads to approximate cancel- 
lations when I attempt to compute reweighted quantities contained in the numerator 
and denominator of ()4.59|) . More details will be given below. 

As we saw in Sec. I4.4.4| det M vanishes for the 5*0 = configurations that were 
discussed above. The phase- quenched quantity {(p^)pg takes this suppression into 
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Figure 3: Phase-quenched averages (hues to guide the eye). 

account. I therefore beheve it gives some hints as to how the CKKU proposal for a 
dynamical lattice spacing works out. In Sec. 14.5.41 I return to the matter of phase 
reweighting. 

4.5.3 Results for the phase-quenched ensemble 

In Fig. 13 (0m'') is obtained in the phase-quenched ensemble; that is, the fermion 
determinant was taken into account in Metropolis updates, using the linear variation 

6 In I det M| = Re Tr [M-^SM] (4.60) 

I remind the reader that larger values of = 1/ go? correspond to finer lattices, 
and thus extrapolate to the continuum limit. However, to hold L = Na fixed, succes- 
sive trajectories must be compared as versus 2N, since is doubled. It can be 
seen from Fig. 01 that for < 6 it is not possible to reliably extrapolate to the large 
L behavior. However, the ga^ = 1/16 results do look encouraging. 

In Fig. m I provide a comparison between the phase- quenched and fully-quenched 
results. It can be seen that the effect of the fermion determinant is dramatic. For 
the finest lattice that I have available, g~^ = 16, the phase- quenched results are 
more than an order of magnitude larger when we look at = 6. Furthermore, the 



51 



A 



-e- 

V 



1.00 



0.50 



0.20 



0.10 



0.05 



0.02 




Figure 4: A comparison of some of the averages obtained in phase-quenched (sohd 
hues to guide the eye) versus fully-quenched (dashed lines to guide the eye) simula- 
tions. Note that the vertical axis in on a log scale. 

phase-quenched results are always of the same order of magnitude as the classical 
estimate ^i^. 

4.5.4 Phase reweighting 

As was reported in [105], and mentioned above, the distribution of a = argdetM 
is essentially flat with respect to the phase- quenched ensemble. In Fig. |21 1 present 
results for a 2 X 2 lattice. It can be seen that it will be quite difficult, if not impossible, 
to obtain a reliable estimate of {exp^ia))^^ . 

E.g., I find that for N = 2, = 2, with 1000 samples, the "noise-to-signal" ratio 

is 

c^(cosa)^ , /(cosQ;)pg = a(sina>^ ^ /(sinQ;)p^ = C(l). (4.61) 

For = 4, = 2, with 1000 samples, I find 

c^(cosa)^ , /(cosa)pg = (T(sina)^ ^ /(sina)pg = 0(10). (4.62) 

On general grounds, these noise-to-signal ratios are expected to get exponentially 
worse as A^ is increased. Similarly discouraging results are found for the other quan- 
tities that need to be estimated in order to make use of the reweighting identity ()4.59|1 . 
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Since reweighting does not even work for small N, it looks like a futile approach by 
which to study the full theory. 

However, the complex phase is not present in the target continuum theory; so, 
it should be possible to address the phase problem. As described above in ^3.8| 
simply taking a — > is not sufficient at a nonperturbative level. The irrelevant 
operators that violate the self-conjugacy of the fermion action are no longer suppressed 
for modes with momentum p ~ 1/a. For this reason, the phase of the fermion 
determinant will always fluctuate wildly. I am currently studying ways to improve 
on the CKKU construction and eliminate the complex phase problem based on this 
new understanding. 

5 Conclusions 

As was explained in ^l.ll many motivations for the study of supersymmetric field 
theories by lattice techniques exist. In ^ we saw that many approaches to lattice 
supersymmetry have been attempted. The plethora of formulations is due to the fact 
that no one of them represents an entirely satisfactory general method. 

In the remainder of this report I have reviewed one of the more promising ap- 
proaches, the one based on dimensional deconstruction. Two 2d SYM theories were 
described, as examples. The (2,2) model described in ^exhibits the basic CKKU 
technique for construction of the class of models. I reviewed the complex phase diffi- 
culty that occurs for the fermion determinant. This problem also arised in the (4,4) 
model of ^ I described how this phase originates from irrelevant, 0{pa) suppressed 
operators that are not present in the continuum. I further explained that since the 
fermion determinant includes modes up to p = 0(a~^), the complex phase does not 
disappear in the a — limit. 

For the (4,4) model that was described in m I studied the classical moduh space 
of the undeformed theory. I then examined the properties of the fermion determinant 
on this moduli space. It was found that this determinant suppresses one branch 
of the moduli space — the one that would have vanishing effective lattice spacing — 
but does not suppress another branch — one that would have nonvanishing effective 
lattice spacing. I showed that the class of minima of the deformed bosonic action is 
gauge-equivalent to the a-configuration, modulo boundary effects. 

Attempts to study the effective lattice spacing by Monte Carlo simulations were 
briefly described. It was seen that the phase-quenched averages — which include the 
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suppression due to the fermion determinant — support the basic picture of CKKU: 
an effective lattice spacing survives when the deformation is scaled to zero in the 
thermodynamic limit. Reweighting by the phase proved to be hopeless, due to a 
nearly fiat distribution. 

It is my hope to present further results regarding renormalization of the effective 
lattice space, using the semiclassical approach, in a future publication. I am cur- 
rently studying formulations that will overcome the phase problem, based on the new 
understanding of its origin that was presented in ^3.81 

As mentioned in ^2.61 I am currently investigating the Sugino construction of 
lattice SYM [52-55]. In a forthcoming article I will give details regarding the lack 
of self-conjugacy for the fermion action, due to irrelevant operators, quite similar to 
what occurs for the CKKU models. I also hope to show how this problem can be 
alleviated through improvements to the action. 

Lastly, I mentioned in §2.91 that I am currently investigating the relationship 
between nonlocal invariances of the lattice partition function and symmetries of the 
continuum theory. Here the issue is that an invariance that touches all time slices 
simultaneously does not have a simple interpretation as the action of an operator that 
commutes with the transfer matrix, and hence hamiltonian. What is needed is to show 
that the operator associated with the transformation localizes in the continuum limit. 
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A A toy (2,2) theory 



In this section I clarify aspects of the "(2,2)" supersymmetry notation, using the 
minimal (2,2) free scalar multiplet theory. There, the Euclidean action takes the 
form 

S = j (fx [-A^dAcp -FF + 2#+a,V'+ - 2#_a^7/^_] . (A.l) 

Here the subscripts on the fermions correspond to an axial U{1) that leaves the scalars 
neutral: 

^±^e±^>±, ^±^e^^>±. (A.2) 

With periodic boundary conditions for all fields, the action is invariant under the 
infinitesmal supersymmetry transformations 

50 = e~ij_ + e+ilj+, SiIj+ = -2ie^d^(j) + e~F, 
5i)_ = 2ie-(9^0-e+F, SF = 2ie~d^ip+ + 2ie^c^'4j_, 

S^p. = -2ie'd,<j)-t^F, SF = 2ie-d,ip+ + 2ie^djip-. (A.3) 
On the other hand, if we make the e's position dependent, 

^S = ^J d^x [d^e-S- + S-cke- + d,e+S+ + S+d,e+] , (A.4) 

where 

5_=z^_<9,0, S^=^.d,<f>, 5+ = V+^0, S+=^+d^(P. (A.5) 

It can be seen from ()A.2|) that S^, and hence the corresponding superalgebra, is of 
left chirality, and that S+ is of right chirality. 
Also note that the equations of motion 

9^^+ = d^^p+ = d^ijj. = djijj. = d^dj(j) = d^d^(f) = (A. 6) 

imply 

= ^+(^), V^- = ^-(^), = 0+(^) + <P~{z), 

tlj+ = i>+{z), tp.=i>.{z), = 0+(z) + 0_(2;). (A.7) 
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Consequently the supercurrents become 



S+iz) = S+{z) = M-^)d^M-^)- (A.8) 

Not only is of left chirality, it is left-moving, since it depends only on z. Similarly 
5+ corresponds to a right-moving superalgebra, depending only on z. Thus, the (2,2) 
notation corresponds equivalently to counting the number of left- and right-moving 
supercurrents. Similar remarks apply for more general values (m,n). 

B Boundedness of the partition function 

One might worry that due to the noncompact classical moduli space of the undeformed 
theory, its partition function (for the purely bosonic theory) is not well-defined. Cer- 
tainly this is the case for case for the Od reduction of d < 4 SU{2) pure Yang-Mills; 
however, this is not the case for the Od reduction of d > 5 SU{2) pure Yang-Mills; 
see Eq. (24) of [101]. This is because, properly speaking, the classical moduli space 
is a set of measure zero in the field integration of the partition function. For li > 5, 
"entropic effects of the measure... overwhelm the possible divergences" [101]. 

I am not aware of any exact result showing that the partition function of the 
CKKU quiver system obtained from the 6d Od SU{kN^) pure Yang- Mills is di- 
vergent.^^ In my view it is unlikely, just from the result that 6d Od SU{2) pure 
Yang-Mills has a finite partition function. Furthermore, if there was a problem with 
the partition function of the undeformed theory, I would have expected an uncon- 
trolled dispersion in the observable ()4.56|) when I attempted to measure it in the 
undeformed theory. This is because flat directions not suppressed by entropic effects 
would allow for the Monte Carlo simulation to wander wildly throughout configura- 
tion space, leading to results that would not converge to reliable average values. That 
this did not occur is further evidence that the partition function is well-defined. Fi- 
nally, the deformation Ssb can be regarded as a regulator of any possible divergence 
associated with the noncompact fiat directions. I obtained stable, identical results as 
the deformation was removed. This provides further evidence that the expectation 
values for the undeformed theory are reliable, and that the quantum moduli space is 
compact. 

^^The diagonal C/(l) from U{kN^) = U{1) x SU{kN^) decouples; its (divergent) contribution to 
the partition function trivially factors out. 



56 



C Deformed fermion matrix 



Suppose we deform the fermion matrix M according to = M + el^^^, where (Im 
is the order (dimension) of the fermion matrix M. Then if Mip = Xip determines an 
eigenvalue A of M, it follows that M^ip = X^ip with = A + e. That is, the set of 
eigenvectors of is identical to that of M ; and, the set of eigenvalues of is just 
the set of eigenvalues of M uniformly shifted by e. 

I now consider e ^ 0^. Suppose that = dimkerM, the number of linearly 
independent eigenvectors of M with eigenvalues of zero. Let det M = Ylx^o denote 
the product of nonzero eigenvalues of M. In the e — 0^ limit, I note 



det M, = J] A, = e"^ J] A (l + ^) = e"^ det M M + ^ ^ + 0{e/Xf j 



• (c.i) 



At^O \ A^O 

Thus there exists a number a = 1 + 0(e/|A|min) such that 

|detM,-e"^detM| < a{dM - n,){e/\X\min), (C.2) 

where \X\min is the minimum nonzero eigenvalue magnitude; i.e., |A|mm = infA^o \M- 
It follows that 

lim — det M, = det M. (C.3) 

In practice I take to be the number of ever-present zeromodes. Then through 
flC.3|) I obtain the product of eigenvalues in the subspace with these modes projected 
out. This is equivalent to having started with an SU{kN'^) mother theory; for exam- 
ple, (TT^ in the (2,2) 2d theory. 

D Conjugation of the action and properties of the 
fermion determinant 

D.l An example from kindergarten 

To begin, I discuss the simple example of single-component Grassmann variables ip, ip, 
with no site indices (i.e., a single-site lattice). Let the partition function have the 
simple form: 

Z{M) = I dijdilje-'^^^ , (D.l) 
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where M is just a complex number. Then according to the usual rules of Grassmann 
integration, Z{M) = M. Now define hermitian conjugation to reverse the order of 
Grassmann variables without the introduction of a sign, so that 

{Z{M)y = iZ{M)y = J d^U^^e-^'^^*^\ (D.2) 

Note that the Grassmann variables are just dummy variables of integration, so that 
we are free to replace 

ip^ ^ip, ^ ip. (D.3) 

Thus we find {Z{M))* = Z{M*) = M*, which is just what one should expect. That 
is, the rules of hermitian conjugation outlined above are self-consistent. It may seem 
strange to be so very explicit about this conjugation rule. However, I could take 
any other redefinition of the integration variables and obtain the same result. For 
example, instead of ()D.3|) I could take 

Then using the Grassmann rules d{uijj) = {l/u)dtjj, etc., 

(Z(M)y = — f ##6"™'^*^*'^ = —uvM* = M*. (D.5) 
uv J uv 

The rule ()D.3|) is just a convenient choice because it restores the measure and the 
action to its original form, apart from the factor M, which has been replaced by M*. 
A more succinct description of the generalized conjugation is as follows. We define 
an operator G that is the composition of hermitian conjugation and ()D.4|) . That is, 

e(v^) = uij, e(^) = vij. (d.6) 

G also has the property that it conjugates complex numbers, transposes matrix ex- 
pressions and reverses the order of fermions. In more sophisticated examples below, 
this freedom of choosing a 6 that defines a generalized conjugation rule will become 
important. In practice B that is unitary and satisfies 0^ = 1 is the most useful. 

Everything that follows will be a generalization of this procedure, done in such a 
way as to remain self-consistent in the above sense. Note that if M* = M, then Z is 
also real. Thus, a hermitian action, 

{i/jMipy = ^^M*i>^ ^pMijj, (D.7) 

using ()D.3|) and M* = M, implies a real partition function Z. 
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D.2 n-component complex Grassmann 

Here I will gather site indices and any spin or internal indices into a collective index 
i = 1, . . . ,n. The partition function is just 

■ ■ ■ #„c^^ne"^*^^, (D.8) 

where M is an n x matrix. As usual, Z{M) = detM. Generalizing ()D.3jl to 

ijl^ij^, i^l^^i, (D.9) 

we find of course that 

(Z(M))* = {Z{M)y = Z{M^) = det = (det M)*, (D.IO) 

which is right. Thus if M"^ = M, one finds that Z{M) is real. Using ()D.9|) . this is 
equivalent to the statement that if the action is self-conjugate, Z{M) is real. 

This self-conjugacy is sufficient but not necessary, as I now show. Suppose that 
instead of (ID. 91) I take 



i)U, vij, (D.ii) 

with U, V SU (n) matrices. Then 

^ (det U det V)-^ j ■ ■ ■ #„#„e-'^^^^^'^. (D.12) 
Because f/, V are unimodular, the prefactor disappears. Thus we find 

det M* = Z(M)* = det(f/MV). (D.13) 
Suppose that we are able to find f/, V such that 

f/MV = e^'^M. (D.14) 

Then it follows that 

77/ 

argdetM = if mod vr. (D.15) 

2 
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As in the 1-component case that was discussed in ^D.ll there exists a more concise 
formulation of this result in terms of a generalized conjugation operator G. That is, 
define 



(D.16) 



Then taking into account ()D.14|) . 



(D.17) 



and ()D.15|) follows. In particular, if the fermion action is self-conjugate in the gener- 
alized sense that 



it follows that det M is real. 

D.3 Imaginary- time continuation of hermit icity 

The most common form of G is the one that follows from hermiticity of the real-time 
action. To extract the imaginary-time continuation of this conjugation operator, 
one just takes Ie = —Hm, etc., in the real-time hermitian conjugation rules. For 
example, since hermitian conjugation treats Minkowski time t^/ as a real variable, 
it transforms Euclidean time as t^; —tE- An example is given in ^3.81 for the 
(2,2) model; in particular, see ()3.57|) . If the real-time action is hermitian, then the 
imaginary-time action is self-conjugate w.r.t. this O. The fermion determinant is real 
in either description of the Grassmann variables. 

E Af = 4: moduli space 

Here I establish the well-known solution to ()4.13|) . One way to see this is as follows 
[109]. First 1 note that 5*0 reduced to the zero modes, which I write as 5*^, takes the 
form 



(D.18) 




(E.l) 
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Now note that the U{1) parts of x, y, z do not appear and can take any value. Thus 
I can restrict our attention the the SU{k) parts, which I choose to express in terms 
of Hermitian matrices ap,bp, p = 1, 2, 3: 



yCJC 



{al + ibl)T'' = ai + ibi, 
{al + ih''^)T^ = a2 + ib2, 
(4 + ib''^)T'' = as + ibs. 



(E.2) 



Substitution into ()E.1|) and a bit of algebra yields 

2 xit^p'^p] +^(K^^i] + i^p^(^<i]y 



9^ 



p,<} 



+ ^{[ap,ag] - [bp,bg]) 
p,g 

i2 



9 



J2 Tr {[ap, a/ + [bp, b,]^ + 2[ap, b,]') 



(E.3) 



p,q 



Using positivity arguments quite similar to those above, one finds that 5*^ > and 
that = iff 



„ ag] = [bp, bg] = [ap, bg] = 0, Vp, q. 



(E.4) 



which is nothing other than Eq. (53) of [109]. Since the matrices are all hermitian 
and they all commute, it is obviously possible to choose a basis which simultane- 
ously diagonalizes them. This basis will be related to the one used in flE.2|) according 
to T'^ — > T'^ = aT^a^ , which is nothing other than the global gauge transforma- 
tions duni). 
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